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1 .  INTRODUCTION 


The  overall  goal  for  this  project  has  been  an  experimental  desiyn  for  the 
determination  of  the  complex  indices  of  refraction  for  small  particles  over  a 
wide  range  in  temperature  and  wavelength.  Two  temperature  ranges  of  interest 
exist:  300  to  500  K  and  1000  to  3000  K.  The  wavelength  range  is  0.3  to 

10.0  pm,  and  the  species  of  interest  are  AI2O3,  Zr02,  MgO,  C,  and  ZrC.  Parti¬ 
cle  sizes  will  range  from  0.1  to  10.0  pm  in  radius.  The  area  of  emphasis  for 
the  experiment  is  the  higher  temperature  range  (1000  to  3000  K)  and  the  infra¬ 
red  wavelengths  (1.0  to  10.0  pm). 

The  project  can  be  segmented  into  three  distinct  tasks:  1)  production  of 
the  experimental  environment,  2)  optical  data  analysis  and  conversion  to  com¬ 
plex  indices  of  refraction,  and  3)  the  optical  diagnostic  design.  These  tasks 
are  discussed  in  detail  in  Sections  2,  3,  and  4,  respectively.  Task  1,  envi¬ 
ronment  production,  considers  the  technologies  available  for  producing  parti¬ 
cles  or  jjarticle  dispersions  at  the  required  temperatures  while  maintaining  an 
appropriate  optical  access.  Task  2,  data  analysis  and  conversion,  considers 
in  detail  what  optical  measurements  are  appropriate,  investigates  their  sensi¬ 
tivity  to  the  complex  index  of  refraction,  defines  a  methodology  for  inverting 
experimental  data  and  identifies  the  errors  associated  with  the  experimental 
measurement  and  inversion  process.  B’inally,  Task  3  details  the  specific  opti¬ 
cal  measurement  system  that  will  produce  the  data  required  by  Task  2. 
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2.  EXPERIMENTAL  ENVIRONMENT  PRODUCTION 


The  creation  of  an  experimental  environment  suitable  for  this  project  is 
addressed  by  defining  appropriate  techniques  for  producing  the  high  and  low 
temperature  environments,  and  by  designing  a  system  that  will  produce  particle 
dispersions  within  this  environment.  Section  2.1  discusses  the  high  tempera¬ 
ture  system,  Section  2.2  presents  the  basic  concepts  for  a  low  temperature 
system,  and  Section  2.3  presents  the  design  for  the  particle  delivery  system. 

2.1  High  Temperature  System 

A  major  portion  of  the  experimental  effort  is  the  production  of  a  steady 
and  definable  high  temperature  environment  for  the  particles.  The  temperature 
range  must  be  1000  to  3000  K,  with  good  optical  access  to  the  temperature 
field.  Previous  experimental  efforts  have  utilized  laser  heating,  direct 
flame  contact  heating,  and  shock  tubes.  In  addition,  the  recent  development 
of  graphite  tube  furnaces  makes  the  3000  K  temperature  attainable  in  a  furnace 
system. 

Electrostatically  levitating  a  single  particle  while  irradiating  it  with 
high  power  density  Nd : YAG  lasers  was  explored  by  Dowling  and  Randall 
(Ref.  1).  This  technique  proved  to  be  exceptionally  difficult  to  implement 
due  to  two  significant  problems.  The  laser-particle  interaction  '.isplaced  the 
particle  out  of  the  levitating  electric  field.  In  addition,  optical  emission 
and  extinction  signals  were  small,  which  mandated  data  collection  and  averag¬ 
ing  over  a  long  time  period.  Subsequent  work  by  Pluchino  et  al.  (Ref.  2)  has 
dropped  the  laser  heating  concept  and  produced  room  temperature  refractive 
index  data  for  carbon  at  0.488  ^m. 

Particle  heating  by  direct  contact  with  a  flame  has  been  utilized  by 
Mularz  and  Yuen  (Ref.  3)  and  Adams  (Ref.  4).  Both  works  were  concerned  with 
temperatures  between  2000  and  3000  K  and  used  alumina  as  the  subject  specie. 
Several  problems  are  inherent  to  this  technique.  The  interaction  of  the 
particle  with  active  species  through  the  flame  zone  will  possibly  produce  a 


contaminated  particle,  thus  adding  a  very  liard  to  quantify  error  to  the  exper¬ 
iment.  In  addition,  the  temperature  of  the  particles  will  change  with  their 
position  in  the  flame,  which  implies  that  a  temperature  map  for  the  flame 
would  have  to  be  produced  and  used  to  determine  the  particle  temperatures. 
Mularz  and  Yuen  (Ref.  3)  also  report  difficulties  associated  with  the  mass 
loading  of  the  flame  with  the  subject  specie  as  well  as  agglomeration  within 
and  downstream  of  the  flame  zone. 

Carlson  (Ref.  5)  also  investigated  a  direct  flame  contact  approach  with  a 
H2-O2  rocket  motor  seeded  with  alumina  and  magnesium  oxide.  Although  he  was 
successful  in  quantifying  the  emissivity  of  AI2O3  as  it  changes  phase,  speci¬ 
fying  the  temperature  field  was  a  significant  problem,  and  background  emission 
from  the  combusting  species  and  products  significantly  complicated  the  optical 
measurements . 

Graphite  tube  furnaces  provide  a  reasonable  method  of  producing  tempera¬ 
tures  up  to  3000  K.  However,  providing  optical  access  and  a  stable  particu¬ 
late  dispersion  would  be  a  major  technical  challenge.  Optical  access  is  dif¬ 
ficult  due  to  the  cooling  required  to  allow  windows  to  survive  the  extreme 
temperatures.  In  addition,  the  background  radiation  due  to  the  furnace  walls 
would  be  very  large.  In  brief,  a  high  temperature  furnace,  although  appearing 
simple  in  concept,  would  be  exceptionally  hard  to  integrate  into  this 
experiment . 

High  temperature  measurements  are  commonly  performed  in  shock  tubes  for  a 
variety  of  experimental  studies.  Nettleton  (Ref.  6)  has  done  an  extensive 
review  of  previous  shock  tube  work  with  non-gas  phase  species  and  discusses  a 
wealth  of  literature  on  the  shock  heating  of  particles,  blast  waves,  explosive 
sensitivity,  slab  combustion,  and  detonation.  Konopka  et  al.  (Ref.  7),  also 
described  in  Calia  el  al.  (Ref.  8),  applied  shock  tube  technology  to  quantify 
the  radiation  from  typical  rocket  plume  particulates  and  reported  the 
technique  as  a  successful  means  of  producing  a  high  temperature  field  of 
particles. 
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In  summary,  high  temperature  experiments  demand  a  well  quantified  and 
steady  environment,  and  three  candidates  have  been  examined  in  the  literature 
in  addition  to  the  graphite  furnace  technique.  Laser  heating  single  i>articles 
appears  to  be  an  exceptionally  difficult  experiment  which  would  require  a 
significant  development  effort.  Flame-particle  interaction  techniques  produce 
particles  of  questionable  composition,  thus  introducing  an  unacceptable  error. 
The  high  temperature  furnace  concept  is  complicated  by  the  problems  of  optical 
access. 

We  believe  a  reflected  zone  shock-tube  technique  provides  a  simple  means 
of  heating  the  particles  and  is  the  most  appropriate  experimental  technique 
available.  In  addition  to  the  work  of  Konopka  et  al.  (Ref.  7)  in  which  they 
reported  no  difficulties  with  the  operation  of  the  shock  tube,  PSI  has  per¬ 
formed  several  shock- tube  studies  and  is  confident  of  the  technique.  Produc¬ 
tion  of  the  high  temperature  environment  will  be  accomplished  most  effectively 
in  a  shock  tube . 

2.1.1  Shock  Tube  Basics  -  Shock  tubes  are  used  to  provide  a  transient  high 
temperature  environment;  their  major  limitation  is  the  time  duration  of  the 
elevated  temperature  field.  The  time  required  for  this  experiment  is  the  time 
necessary  for  temperature  equilibration  and  subsequent  optical  measurements. 
This  "working"  time  is  dependent  upon  several  of  the  parameters  relevant  to 
the  shock  tube's  operation.  In  this  section  a  development  of  the  specific 
dependencies  is  presented  along  with  a  brief  discussion  of  the  basics  of  shock 
tube  operation. 

Shock  tubes  are  typically  discussed  with  x-t  diagrams,  as  in  Leipman  and 
Roshko  (Ref.  9)  and  Gaydon  and  Hurle  (Ref.  10),  which  allow  description  of  the 
evolution  in  distance  and  time  of  the  shock  wave,  contact  discontinuity,  and 
rarefaction  wave,  as  shown  in  Figure  1.  Before  the  diaphragm  breaks,  at  time 
0,  the  driver  section  is  filled  with  quiescent  gas  at  pressure  P4  and  tempera¬ 
ture  T4.  Similarly  the  test  section  is  filled  with  quiescent  gas  at  pressure 
P)  and  temperature  Tj .  When  the  diaphragm  breaks,  a  shock  wave  propagates 
through  the  test  section  followed  by  the  contact  discontinuity,  which  is  a 
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STATION 

(a) 


Figure  1(a).  Fluid  motion  in  a  shock  tube. 
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surface  that  separates  the  test  gas  and  driver  gas.  Pressure  and  velocit’ 
are  equal  across  this  surface.  A  rarefaction  front,  traveling  at  the  local 
speed  of  sound,  propagates  into  the  driver  section  and  is  followed  by  an 
expansion  fan  which  is  a  field  of  varying  temperature,  pressure,  and  velocity 
that  separates  the  quiescent  gas  in  Region  4  from  the  flowing  constant  temper¬ 
ature  and  pressure  gas  in  Region  3.  Thus,  when  the  diaphragm  breaks,  a  tran¬ 
sient  flow  begins.  Regions  1  and  4  are  the  original,  unaffected  regions  in 
the  shock  tube.  Regions  2  and  3  are  the  flowing  regions,  with  equal  velocity. 
Pressures  and  temperatures  within  the  shock  tube  at  a  fixed  time,  tj,  are 
shown  in  Figure  1(b).  As  previously  stated,  the  pressures  and  temperatures  in 
Regions  1  and  4  are  at  their  original  levels.  Region  2  is  the  gas  that  has 
been  heated  by  the  shock,  explaining  the  increased  pressure  and  temperature. 
Region  3  contains  the  gas  that  has  undergone  an  expansion,  thus  lowering  the 
temperature  and  pressure  from  the  starting  value  in  Region  4. 


The  optical  measurements  for  this  project  will  take  place  in  Region  5 
which  is  a  quiescent  region  through  which  the  reflected  shock  has  passed.  The 
problem  from  an  experimental  design  point  of  view  is  to  predict  the  tempera¬ 
ture,  pressure  and  duration  of  Region  5  to  insure  that  these  parameters  are 
appropriate  for  the  experiment.  The  equations  which  describe  this  behavior 
are  straightforward  algebraic  expressions  and  are  derived  in  Gaydon  and  Hurle 
(Ref.  10).  However,  their  form  dictates  iterative  solutions.  The  two  rela¬ 
tions  required  to  specify  the  pressure  in  the  reflected  shock  zone,  P5,  are 
given  by  Eqs.  (1)  and  (2). 

V1 

4  1  '  ,L  ♦  lvi)<e-i)  J 

(3V1)  p/p,  -  yn 

V  2  ■  (yT+1)  +  (yt-1>  VP1 

liquations  (3),  (4),  and  (5)  describe  the  temperature  in  Region  5. 


(yd-0(V«4)<m> 
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£  =  V?! 


3a 


£  '  -  P5/P2 


f({)  -{ivnr 

T5/T 1  =  f(c)  f(c • ) 


(yt-i)£  +  (yt+i) 


yt+iJ£  +  (yt-1j 


A  full  list  of  variables  and  their  meanings  is  given  in  Table  1 
TABLE  1.  Variables  for  the  Shock  Tube  Equations. 


a 

-  speed  of  sound 

L 

-  length 

M 

-Mach  number 

P 

-  pressure 

u 

-  velocity 

X 

-  distance  from  the  diaphragm  station 

t 

-  time 

T 

-  temperature 

Y 

-  specific  heat  ratio 

Subscrip 

ts 

c  - 

contact  surface 

D  - 

driver  section 

T  - 

test  section 

r  - 

rarefaction  front 

1,2, 3,4  - 

1 _ 

conditions  shown  in  Figure  1 

The  working  time  in  Region  5  is  determined  by  the  interaction  oi  the 
reflected  shock  with  the  contact  surface,  provided  that  the  driver  meets  a 
certain  minimum  length  requirement.  The  shock-contact  surface  interaction  is 
shown  in  Figure  2  along  with  the  working  time,  tw,  which  is  the  time  between 
wave  arrivals  at  the  end  wall. 

As  shown  in  this  figure,  the  reflected  shock  and  contact  surface  meet  and 
produce  two  waves:  a  weak  reflected  shock  and  a  transmitted  shock.  The  work¬ 
ing  time  in  the  reflected  shock  zone  is  the  sum  of  two  intervals:  the  time 
for  the  reflected  shock  to  meet  the  contact  surface,  t ) ,  and  the  time  for  the 
shock  reflected  off  the  contact  surface  to  strike  the  end  wall,  i2«  The  first 
interval  is  found  by  equating  the  position  of  the  reflected  shock, 

Eq.  (6),  with  that  of  the  contact  surface,  Eq.  (7),  to  find  the  time,  t, 

Eq.  (8).  Ti,  Eq.  (10),  is  the  difference  between  this  time  and  the  time  for 
the  incident  shock  to  hit  the  end  wall. 


i 


Figure  2.  Interaction  of  the  reflected  shock  with  the  contact  surface. 
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xc  =  u3t 


(6) 


xr  =  L2  -  (M2a2  -  u2)(ti-t') 


(7) 


t. 

1 


L0  +  t1  [  M 


V2 


(8) 


f  =  l2/m,a, 


(9) 


L_  Ma  „  +  M  a ^  -  u„ 
2  r  1  1  2  2  2 1 

T 1  '  vT  ~ 


M„a„ 
2  2 


(  10) 


The  second  interval,  x2>  is  conservatively  approximated  by  assuming  the 
weak  reflected  shock  to  have  a  Mach  number  of  1.  Therefore,  t2  is  equal  to 
the  distance  traveled  by  the  reflected  shock  divided  by  the  speed  of  sound  in 
the  reflected  shock  zone,  Eq.  (11),  and  the  total  working  time,  xw,  is  the  sum 
of  t i  and  t  2,  Eq.  ( 12) . 


M_a 


_  -  u  ] 

2  2  2J 


(11) 


Tw  =  T  1  +  T2 


(12) 


This  discussion  of  the  working  time  in  the  reflected  shock  zone  assumes 
that  the  reflected  head  of  the  expansion  wave  has  not  passed  through  the  con¬ 
tact  surface  and  is  a  reasonable  assumption  if  the  driver  section  is  long 
enough.  Figure  3  illustrates  the  evolution  of  the  rarefaction  wave  for  the 
minimum  length  driver  section  required  for  no  rarefaction  effects  in  the 
reflected  shock  zone  during  the  working  time,  tw.  This  condition  dictates 
that  the  rarefaction  head  meets  the  contact  surface  no  sooner  than  the  end  of 
the  period  x  i .  The  time  at  which  the  expansion  front  meets  the  contact  sur¬ 
face  is  given  by  Eq.  (13);  the  development  of  this  equation  is  given  in  Gaydon 
and  Hurle  (Ref.  10).  This  time  is  equated  to  the  time  it  takes  the  reflected 
shock  to  meet  the  contact  surface,  Eq.  (8),  which  yields  the  required  ratio  of 
driver  length  to  test  length,  Eq.  (14). 
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Figure  3.  Rarefaction  wave  interaction  with  the  contact  surface. 


t. 
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L  /L 
D 


T 


Y  +1 
1  D 

W) 


M  a 
11 


+  V2  "  U2J 


M  a 
11 


M  a„ 
2  2 


+ 


-( 


1 


(13) 


(14) 


Given  these  relationships,  the  working  time  in  the  shock  tube  can  be 
related  to  the  test  gas  and  the  temperature  in  the  reflected  shock  zone,  as 
tabulated  for  hydrogen-argon,  helium-argon  and  nitrogen-argon  systems  in 
Tables  2,  3  and  4.  These  tables  clearly  indicate  that  the  v/orking  time 
decreases  at  the  higher  temperatures.  The  driver  length- to- test  length  ratio, 
listed  in  the  tables  illustrates  that  the  maximum  driver  length  is  required  at 
the  lower  temperatures.  The  initial  and  reflected  pressure  ratio  in  the  shock 
tube,  P4/P 1 ,  and  P5/P1,  are  also  given  in  Tables  2,  3,  and  4  as  a  function  of 
temperature  in  the  reflected  shock  zone. 
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TABLE  2.  Shock  Conditions  for 


a  Hydrogen  Argon  System. 


Working  time,  tau,  is  given  in  seconds  per  meter  of 
LDLT  is  the  minimum  ratio  of  driver- to- test 
Ml  and  M5  are  the  incident  and  reflected  shock 


test  section 
lengths, 
strengths. 


INITIAL  GAS  TEMPERATURE  IS  JOO. 
ORIVER  GAS  IS  HYDROGEN  .  GAMMA  IS 
TEST  GAS  IS  ARGON  GAMMA  IS 


T3  TAU  LDLT 

BOO.O  0. 131E-02  1.690799 

330.0  0.I22E-O2  1.614811 

’00.0  0.113E-02  1.347609 

*30.0  O.I06E-02  1.487687 

1000.0  0.993E-03  1.433874 

1030.0  0 . 940E-03  l . 383234 

1100.0  0.892E-O3  1.341012 

1130.0  0.848E-O3  1.300604 

1200.0  0. 810E-03  1.263303 

1230.0  0. 773E-03  1.229298 

1300.0  0. 743E-03  1.197636 

1330.0  0. 713E-03  1.168227 

1400.0  0.688E-03  1.140823 

1450.0  0.663E-03  1.113204 

1300.0  0.643E-03  1.091199 

1330.0  0.622E-03  1.068643 

1600.0  0.603E-03  1.04779/ 

1630.0  0.3B6E-03  1.027347 

1700.0  0.370E-03  1.008379 

1730.0  0.333E-03  0.990403 

1800.0  0.340E-03  0.973341 

1830.0  0.327E-03  0.937118 

1900.0  0. 313E-03  0.941666 

1930.0  0.303E-03  0.926924 

2000.0  0. 492E-03  0.912848 

2030.0  0. 481E-03  0.399379 

2100.0  0.471E-03  0.886489 

2130.0  0.462E-03  0.874124 

2200.0  0.453E-03  0.862260 

2230.0  0. 444E-03  0.830838 

2300.0  0.436E-03  0.839893 

2330.0  0.429E-03  0.829336 

2400.0  0.421E-03  0.819164 

2430.0  0.414E-O3  0.809332 

2300.0  0.407E-03  0.799882 

2330.0  0.401E-03  0.790731 

2600.0  0.393E-03  0.781889 

2630.0  0.389E-O3  0.773330 

2700.0  0.383E-03  0.763047 

2730.0  0.378E-03  0.737021 

2800.0  0.372E-03  0.749241 

2830.0  0.367E-03  0.741694 

2900.0  0.362E-03  0.734368 

2930.0  0. 338E-03  0.727236 

3000.0  0. 333E-03  0.720343 

3030.0  0.349E-O3  0.713624 

3100.0  0.343E-03  0.707083 

3150.0  0.340E-O3  0.7O0723 

3200.0  0.336E-03  0.694331 


1.41,  R  IS  4122.0 
1.67,  R  IS  208.0 


P2/P1  P4/P1  P3/P1  m 

3.490  4.729  9.283  1.729 

3,737  3.203  10.427  1.790 

4.026  3.692  11.606  1.849 

4.293  6.196  12.816  1.906 

4.363  6.714  14.033  1.962 

♦.336  7.246  13.320  2.017 

3.108  7.792  16.610  2.070 

3.380  8.332  17.921  2.122 

5.632  8.923  19.232  2.172 

3.923  9.312  20.602  2.222 

6.198  10.112  21.968  2.270 

6.471  10.723  23.330  2.318 

6.743  11.332  24.743  2.363 

7.018  11.992  26.133  2.411 

7.292  12.644  27.376  2.436 

7.367  13.310  29.008  2.300 

7.841  13.989  30.432  2.343 

8.113  14.680  31.904  2.386 

8.390  13.3B3  33.367  2.628 

3.663  16.103  34.837  2.670 

8.940  16.834  36.316  2.711 

9.213  17.377  37.802  2.731 

9.490  18.334  39.293  2.790 

9.763  •  19.104  40.793  2.83M 

10.040  19.386  42.300  2.368 

10.313  20.682  43.812  2.906 

10.391  21.491  43.329  2.944 

10.866  22.313  46.832  2.981 

11.142  23.148  48.379  3.018 

11.417  23.996  49.911  3.034 

11.693  24.837  31.447  3.090 

11.968  23.732  32.988  3.123 

12.244  26.620  34.332  3.160 

12.320  27.321  36.080  3.193 

12.793  28.433  37.632  3.229 

13.071  29.364  39.187  3.263 

13.347  30.303  60.743  3.297 

13.623  31.260  62.307  3.330 

13.899  32.228  63.871  3.363 

14.173  33.211  63.439  3.396 

14.431  34.206  67.008  3.428 

14.726  33.216  68.381  3.460 

13.002  36.239  70.136  3.492 

13.278  37.276  71.732  3.323 

13.334  38.328  73.312  3.333 

15.830  39.393  74.893  3.383 

16.107  40.472  76.477  3.616 

16.383  41.363  78.062  3.646 

16.659  42.672  79.630  3.677 


M3 

1.326 
1.333 
1.383 
1.608 
1.631 
1.633 
1.673 
1.692 
1.710 
1.726 
1.742 
1.736 
1.770 
1.783 
1.793 
1.807 
1.818 
1.828 
1.838 
1.848 
1.837 
1.363 
1.874 
1.382 
1.389 
1.396 
1.903 
1.910 
1.916 
1.922 
1.928 
1.934 
1.939 
1.943 
1.930 
1.933 
1.939 
1.964 
1.968 
1.973 
1.977 
1.981 
1.983 
1.908 
1.992 
1.996 
1.999 
2.002 
2.006 


TABLE  3. 


Shock  Conditions 


f°r  a  Helium  Argon  System. 


noMT™L  MS  IS  300  i( 

test  mWl'T*  aiiTS  15 ,  J'tis  15  <073-> 

*  ■  o/ >  R  IS  208.0 


T3 

800.0 
830.0 
900.0 
930.0 
1000.0 
1050.0 
1100.0 
1130.0 
1200.0 
1230.0 
1300.0 
1330.0 
1400.0 
1430.0 
1300.0 
1330.0 
1600.0 
1630.0 
1700.0 
1730.0 
1800.0 
1830.0 
1900. 0 
1930.0 
2000.0 
2030.0 
2100.0 
2130.0 
2200.0 
2230.0 
2300.0 
2330.0 
2400.0 
2430.0 
2300.0 
2330.0 
2600.0 
2630.0 
2700.0 
2730.0 
2800. 0 
2830.0 
2900. 0 
2930.0 
3000. 0 
3030.0 
3100.0 
3130.0 
3200.0 


TflU 

0. 131E-02 
0- 1226-02 
0.113E-02 
0. 106E-02 
0. 993E-03 
0. 94OE-03 
0.892E-03 
0.B48E-03 
0.810E-03 
0. 773E-03 
0. 743E-03 
0.713E-03 
0. 688E-03 
0.663E-03 
0.643E-03 
0. 622E-03 
0.603E-03 
0.386E-03 
0. 370E-03 
0.333E-03 
0.340E-03 
0. 327E-03 
0.313E-03 
0. 303E-03 
0.492E-03 
0.481E-03 
0.471E-03 
0. 462E-03 
0. 433E-03 
0.444E-03 
0.436E-03 
0. 429E-03 
0.421E-03 
0.414E-O3 
0.407E-03 
0.401E-O3 
0.393E-03 
0.389E-03 
0.383E-03 
0. 378E-03 
0.372E-03 
0. 367E-03 
0. 362E-03 
0.338E-03 
0.333E-03 
0.347E-O3 
0.343E-03 
0.3406-03 
0. 336E-03 


UJLT 
1.301674 
1.243173 
1.191437 
1.143306 
1. 103878 
1.066432 
1.032387 
1.001279 
0.972718 
0.944383 
0.922008 
0.899367 
0.878270 
0.838348 
0.840067 
0.822702 
0.806344 
0.790910 
0.774307 
0.762470 
0.749333 
0.734844 
0.724948 
0.713399 
0.702762 
0.692393 
0.682469 
0,672931 
0.663817 
0.633039 
0.646397 
0.638470 
0.630639 
0.623003 
0.613794 
0.608730 
0.601943 
0.393333 
0.388977 
0.382798 
0. 376808 
0.321407 
0.318232 
0.313130 
0.312133 
0. 309243 
0. 306409 
0.303633 
0. 300969 


P  2/Pi 
3.490 
3.  737 
4.026 
4.293 
4.363 
4.836 
3.108 
3. 380 
3.632 
3. 923 
6.198 
6.471 

6.743 
7.018 
7.292 
7.367 
7.841 
8.113 
8.390 
8.663 
8.940 
9.213 
9.490 
9.763 
10.040 
10.313 
10.391 
10.866 
11.142 
11.417 
11.693 
11.968 
12.244 
12.320 
12.793 
13.071 
13.347 
13.623 
13.899 
14.173 
14.431 
14.726 
13.002 
13.278 
13.334 
13.830 
16.107 
16.383 
16.639 


P4/P1 
3.619 
6.263 
6.943 
7.632 
8.394 
9.167 
9.972 
10.810 
11.681 
12.384 
13.322 
14.493 
13.498 
16.538 
17.614 
18.723 
19.872 
21.036 
22.277 
23. 337 
24.834 
26.171 
27.347 
28.964 
30.421 
31.921 
33.462 
33.046 
36.674 
38.347 
*0.064 
41.827 
43.637 
43.493 
47.400 
49.333 
31.339 
33.414 
33.320 
37.679 
39.891 
62.137 
64.479 
66.833 
69.289 
71.781 
74.733 
76.943 
79.613 


P3/P1 
9.283 
10.427 
11.606 
12.816 
14.033 
13.320 
16.610 
17.921 
19.232 
20.602 
21.968 
23.330 
24.743 
26.133 
27.376 
29.008 
30.432 
31.904 
33.367 
34.837 
36.316 
37.802 
39.293 
40.793 
42.300 
43.812 
43.329 
46.832 
48.379 
*9.911 
31.447 
32.988 
34.332 
36. 080 
57.632 
39.187 
60. 745 
62.307 
63.871 
63.  439 
67. 008 
68.381 
70.136 
71.732 
73.312 
74.893 
76. 477 
78.062 
79. 630 
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1.729 
1.790 
1.849 
1.906 
1.962 
2.017 
2.070 
2.122 
2.172 
2.222 
2.270 
2.318 
2.363 

2.411 
2.436 
2.300 
2.343 
2.386 
2.62B 
2.670 
2.711 
2.731 
2.790 
2.830 
2.868 
2.906 
2.944 
2.981 
3.018 
3.034 
3.090 
3.123 
3.160 
3.193 
3.229 
3.263 
3.297 
3.330 

3.363  | 

3.396  i 

3.428  l 

3. 460  1 

3. 492  | 

3.523  1 

3.333  i 

3.383  i, 

3.616  |. 

3. 646 
3. 


2. Om 
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2.006 


TALil.K  4 


Shock  Conditions  for  a  Nitrogen  Argon  System 


INITIAL 

GAS  TEMPERATURE 

o 

o 

n 

tn 

K 

DRIVER  GAS  IS  NITROGEN 

.  GAMMA  IS 

1.40. 

R  IS  296 

.6 

TEST  GAS 

IS  ARGON 

GAMMA  IS 

1.67.  R 

IS  208.0 

T3 

TAU 

LDLT 

P2/PI 

P4/P1 

P3/P1 

Ml 

M3 

soo.o 

0. 1 3 IE— 02 

0.270073 

3.490 

11.603 

9.283 

1.729 

1.326 

850.0 

0. 122E-02 

0.247932 

3.737 

13.697 

10.427 

1.790 

1.333 

900. 0 

0. 1 13E-02 

0.228723 

4.026 

16.043 

11.606 

1.849 

1.383 

950.0 

0. 106E-02 

0.211839 

4.293 

18.667 

12.816 

1.906 

1.608 

1000.0 

0.996E-03 

0. 196936 

4.563 

21.583 

14.033 

1.962 

1.631 

1050.0 

0.940E-03 

0. 183693 

4.836 

24.822 

15.320 

2.017 

1.653 

1100.0 

0. B92E-03 

0.171822 

3.108 

28.403 

16.610 

2.070 

1.673 

1150.0 

0. 848E-03 

0.161133 

5.380 

32.339 

17.921 

2.122 

1.692 

1200.0 

0.810E-03 

0.151462 

3.652 

36.715 

19.252 

2.172 

1.710 

1250.0 

0.775E-03 

0. 142672 

3.923 

41.503 

20.602 

2.222 

1.726 

1300.0 

0. 743E-03 

0. 134650 

6.198 

46.737 

21.968 

2.270 

1.742 

1250.0 

0. 713E-03 

0. 127303 

6.4M 

52.312 

23.330 

2.318 

1.756 

1400.0 

0.688E-03 

0. 120550 

6.745 

58.B06 

24.743 

2.365 

1.770 

1450.0 

0.665E-03 

0.114322 

7.018 

65.682 

26. 133 

2.411 

1.783 

1500.0 

0. 643E-03 

0.108564 

7.292 

73. 180 

27.576 

2.456 

1.795 

1550.0 

0. 622E-03 

0. 103225 

7.567 

81.350 

29.008 

2.500 

1.807 

1600.0 

0.603E-03 

0.098263 

7.841 

90.241 

30.452 

2.343 

1.B1B 

1650.0 

0.386E-03 

0.093640 

8.113 

99.905 

31.904 

2.386 

1.828 

1700.0 

0.370E-07 

0.089324 

8.390 

110.402 

33.367 

2.628 

1.838 

1750. C 

0. 555E-03 

0.083287 

8.665 

121.791 

34.837 

2.670 

1.848 

1800.0 

0. 340E-03 

0.081302 

8.940 

134.141 

36.316 

2.711 

1.837 

1850.0 

0.327E-03 

0.077950 

9.215 

147.316 

37.802 

2.751 

1.B65 

1900.0 

0.315E-03 

0.074610 

9.490 

161.995 

39.295 

2.790 

1.874 

1950.0 

0.503E-03 

0.071463 

9.765 

177.660 

40.793 

2.830 

1.882 

2000.0 

0. 492E-03 

0.068496 

10.040 

194.591 

42.300 

2.868 

1.889 

2050.0 

0.481E-03 

0.063693 

10.315 

212. B90 

43.812 

2.906 

1.896 

2100.0 

0. 471E-03 

0.063043 

10.391 

232.638 

45.329 

2.944 

1.903 

2130.0 

0. 462E-03 

0.060334 

10.866 

253. 939 

46.852 

2.981 

1.910 

2200.0 

0. 433E-03 

0.058133 

11.142 

276. 932 

48.379 

3.018 

1.916 

2230.0 

0.444E-03 

0.033898 

11.417 

301.745 

49.911 

3.054 

1.922 

2300.0 

0. 436E-03 

0.033733 

11.693 

328. 458 

31.447 

3.090 

1.928 

2330.0 

0.429E-03 

0.031716 

11.968 

357.239 

52.988 

3.123 

1.934 

2400.0 

0.421E-03 

0.049776 

12.244 

388.224 

54.332 

3.160 

1.939 

2430.0 

0. 414E-03 

0.04792B 

12.320 

421. 3B2 

56.0B0 

3.195 

1.943 

2300.0 

0.408E-03 

0.046167 

12.793 

457.471 

57.632 

3.229 

1.930 

2330.0 

0.  401E-03 

0.044486 

13.071 

496.083 

39. 187 

3.263 

1.933 

2600.0 

0. 393E-03 

0.042BB1 

13.347 

537.583 

60.745 

3.297 

1.959 

2650.0 

0. 389E-03 

0.041347 

13.623 

582.223 

62.307 

3.330 

1.964 

2700.0 

0. 383E-03 

0.039882 

13.899 

630. 176 

63.871 

3.363 

1.968 

2750.0 

0.378E-03 

0.038479 

14.173 

681.711 

63.439 

3.396 

1.973 

2800.0 

0.372E-03 

0.037136 

14.431 

737. 066 

67.008 

3.428 

1.977 

2830.0 

0.367E-03 

0.035849 

14.726 

796.508 

6B.S81 

3.460 

1.981 

2900. 0 

0. 363E-03 

0.034616 

15.002 

860.352 

70.156 

3.492 

1.985 

2930.0 

0.358E-03 

0.033433 

13.279 

928.862 

71.732 

3.323 

1.988 

3000.0 

0. 353E-03 

0.032298 

15.334 

1002.415 

73.312 

3.555 

1.992 

3030.0 

0. 349E-03 

0.031209 

15. 830 

1081.347 

74.893 

3.583 

1.996 

3100.0 

0.343E-03 

0.030162 

16.107 

1166.067 

76.477 

3.616 

1.999 

3130.0 

0.340E-03 

0.029136 

16.383 

1256.924 

78.062 

3.646 

2.002 

3200.0 

0.336E-03 

0.028189 

16.659 

1334.419 

79.630 

3.677 

2.006 

The  relationships  used  to  assemble  Tables  2,  3,  and  4  provide  the  basic 
tools  for  designing  the  shock  tube  experiment.  The  next  step  is  to  define  the 
working  time  required  in  the  reflected  shock  zone. 

2.1.2  Particle-Temperature  Equilibrium  -  Two  parameters  of  importance  when 
considering  the  energy  exchange  between  the  particles  and  the  heated  gas  in 
the  reflected  zone  of  the  shock  tube  are  the  change  in  the  gas  temperature  and 
the  time  it  takes  for  particle-gas  equilibration  occur.  The  change  in  gas 
temperature  is  directly  related  to  the  mass  loading  of  particles  and  the  ini¬ 
tial  temperature  difference  between  the  particles  and  the  gas.  An  upper  bound 
on  number  densities  for  the  experiment  is  provided  by  setting  the  optical 
depth  equal  to  0.1.  The  optical  depth  is  defined  by  Bouguer's  law,  Eq.  (15a), 
and  is  a  means  of  quantifying  the  attenuation  of  energy  in  an  absorbing  and 
scattering  media.  Optical  depth  also  quantifies  the  type  of  scattering  for 
the  media:  optical  depths  less  than  0.1  produce  single  scattering,  optical 
depths  between  0.1  and  0.3  produce  double  scattering,  and  optical  depths 
greater  than  0.3  produce  multiple  scattering.  This  upper  bound  on  number 
densities  is  expressed  in  Eq.  (15b)  which  follows  directly  from  setting  the 
optical  depth  equal  to  0.1  in  the  relationship  between  optical  depth,  particle 
number  density,  extinction  efficiency,  particle  radius,  and  the  optical  path 
length. 

I/IQ  =  e“T  \ 15a) 

0.1 

N  «=  - - -  (15b) 

P  2 

v  tt  l  r  Q  . 

ext 

where : 

I/I0  -  ratio  of  attenuated  to  incident  intensity 
t  -  optical  depth 

i  -  optical  path  length  in  the  attenuating  media 
r  -  partical  radius 

-  extinction  efficiency 
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An  estimate  of  the  number  density  and  the  resulting  mass  loading  can  be  made 
by  using  1.0  as  an  estimate  for  the  extinction  efficiency  (Qext  can  be  as 
great  as  4;  see  Figure  42) ,  relating  tne  resultant  number  density  to  a  mass 
loading  of  the  particles.  This  yields  number  densities  of  3  x  103  to  3  x  107 
particles/cm3,  corresponding  to  the  radii  equal  to  10.0  and  0.1  um»  respec¬ 
tively,  which  imply  AI2O3  mass  loadings  of  5  x  10“-*  to  5  x  10“?  gm/crn3.  This 

mass  of  particles  will  depress  the  initial  temperature  in  the  reflected  shock 

zone.  A  simple  energy  balance  can  be  used  to  quantify  this  temperature 

change,  as  shown  in  Eq.  (16). 


m  c 


where : 

ATg/ATp  -  ratio  of  the  changes  in  the  gas  temperature  to  the  change  in 
the  particle  temperature 

Cp/cg  -  ratio  of  the  particle  specific  heat  to  the  gaseous  specific 
heat 

mp/mg  -  ratio  of  the  particle  mass  to  the  mass  of  the  gas. 

As  an  example,  for  extinction  efficiencies  equal  to  one  with  an  argon  test 
gas,  a  working  pressure  of  20  atm,  and  a  test  temperature  of  2000  K,  the 
change  in  gas  temperature  will  be  between  0.025  and  2.5  percent. 

The  time  required  to  reach  temperature  equilibrium  is  a  more  immediate 
design  factor  for  the  shock  tube  experiment  and  has  been  investigated  in 
greater  detail.  A  lumped  mass  analysis  will  be  used  to  model  the  particle 
temperature  response;  any  temperature  gradients  within  the  particle  will  be 
neglected.  This  approximation  is  commonly  used  for  systems  in  which  the  rate 
of  heat  transfer  in  one  part  of  the  system  is  much  greater  than  that  in 
another.  For  no  convective  heat  transfer  between  the  gas  and  the  particle, 
the  ratio  of  the  gas  to  particle  conductivity,  kg/kp,  is  a  reasonable  measure 
of  the  accuracy  of  this  statement.  In  all  cases  this  conductivity  ratio  will 
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be  less  than  0.1,  thus  justifying  the  assumption.  An  energy  balance  for  the 
system  is  given  by  Eq.  (17).  A  simplifying  assumption,  used  for  direct  solu¬ 
tion,  is  to  consider  the  overall  system  to  be  one  particle  surrounded  by  a 
fluid  at  the  final  equilibrium  temperature,  Te.  This  assumption  allows  direct 
integration  of  Eq.  (17)  by  ignoring  the  variability  in  the  gas  temperature. 

The  error  that  is  introduced  will  be  small  and  will  cause  an  overestimation  of 
the  time  to  reach  equilibrium,  thus  providing  a  slightly  conservative  esti¬ 
mate.  Equation  (18)  is  the  integrated  form  of  the  energy  balance. 

9T 

pcv  ^  «  -  hA( T  -  Tft)  (17) 

"  -***>(-  (,B> 

where : 

0  =  <Tp-Te)/(T0-Te) 

Tp  -  particle  temperature 

Te  -  equilibrium  temperature  of  particle  and  gas 

T0  -  initial  particle  temperature 
h  -  convective  heat  transfer  coefficient 
A  -  particle  surface  area 
p  -  particle  density 
c  -  specific  heat  of  particle 
V  -  particle  volume 

The  integration  for  Eq.  (17)  assumed  a  constant  specific  heat  for  the 
particle.  This  approximation  is  appropriate  provided  the  value  which  is 
actually  used  is  chosen  to  give  a  conservative  value  for  the  equilibration 
time.  A  final  approximation  is  required  for  a  usable  solution;  the  convective 
heat  transfer  coefficient,  h,  is  not  specified  but  may  be  approximated  by 
quick  consideration  of  the  fluid  mechanical  aspects  of  the  problem.  The  rela¬ 
tive  velocity  between  the  gas  and  the  particle  will  be  very  small.  This 
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implies  that  conduction  is  the  primary  mode  of  heat  transport,  which  gives  h 
as  the  ratio  of  gas  conductivity  to  particle  radius.  If  temperature  equilib¬ 
rium  is  defined  as  less  than  two  percent  difference  in  the  particle  and  air 
temperature,  the  equil  ibration  time  is  given  by  Eq.  (19): 


where : 

r  -  particle  radius 

k  -  gas  thermal  conductivity 

Figure  4  displays  this  equilibration  time  versus  particle  radius  at  1000 
and  3000  K  for  AI2O3;  the  smaller  times  for  3000  K  are  due  to  the  increase  in 
thermal  conductivity  of  the  gas.  Large  radius  particles  take  a  significant 
time  to  come  to  temperature  equilibrium,  and  a  means  of  decreasing  this  time 
is  necessary.  A  very  direct  method  of  decreasing  this  time  is  to  change  the 
thermal  conductivity  of  the  gas  by  increasing  the  pressure  in  the  reflected 
zone  of  the  shock  tube.  This  pressure  could  easily  be  driven  up  to  20  atm 
which  would  increase  the  thermal  conductivity  by  a  factor  of  20.  The  result¬ 
ing  equilibration  times  are  also  shown  in  Figure  4,  and  are  within  the  working 
times  for  reasonable  length  shock  tubes. 

Figure  4  also  displays  the  available  time  for  particle  equilibration  in 
the  PSI  4-in.  shock  tube  facility.  The  test  section  is  7.2  meters  long  which 
results  in  the  working  times  shown.  The  driver  section  is  4  meters  in  length 
which  yields  a  driver  to  test  length  ratio  of  0.56.  Tables  2,  3,  and  4  can  be 
used  to  determine  the  appropriate  driving  gas  for  the  temperature  in  the 
reflected  shock  zone.  The  lightest  driving  gas  which  produces  no  unwanted 
rarefaction  effects  in  the  reflected  shock  zone  is  used  whenever  possible  be¬ 
cause  of  the  lower  operating  pressures  associated  with  the  lighter  gas.  Argon 
is  used  as  the  test  section  gas  because  of  its  chemically  and  optically  inert 
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behavior.  Table  5  lists  operating  temperatures  between  1000  to  3000  K  along 
with  the  working  time  available  and  the  shock  tube  gases  to  use. 

2.2  Low  Temperature  Environment  Concept 

A  conceptual  design  was  produced  for  the  low  temperature  measurements  in 
this  study.  This  hardware  can  be  configured  similar  to  a  shock  tube  so  that 
particle  delivery  hardware  and  optical  diagnostics  developed  for  the  shock 
tube  experiment  can  be  readily  adapted  to  the  low  temperature  setup.  A  tubu¬ 
lar  flowing  furnace  system  is  envisioned,  as  shown  in  Figure  5.  The  particle 
seeding  section  is  positioned  upstream  of  the  optical  measurement  section. 
Particles  will  remain  entrained  in  the  flow  due  to  their  small  terminal  veloc¬ 
ities.  These  velocities  may  be  computed  from  a  Stokes  Drag  formulation  and 
the  results  are  presented  in  Figure  6. 

Low  terminal  velocities  are  very  accommodating  since  they  imply  that  once 
established  a  dispersion  of  particles  will  be  relatively  stable.  The  low  tem¬ 
perature  experiment  will  be  temperature-limited  by  the  infrared  background 
emission  attributable  to  the  furnace  system  and  by  material  interaction  prob¬ 
lems  between  furnace  and  the  optical  access  windows.  Tnis  temperature  limit 
will  probably  be  between  600  and  800  K. 

TABLE  5.  Shock  Tube  Summary. 

4-in.  I.D.  Tube 
Test  Section  -  7.2m 
Driver  Section  -  4m 


Working  Time 
(ms) 

Gases-Driver/ 

Test 

1000 

7.2 

Ni trogen/Argon 

1500 

4.6 

Nitrogen/Argon 

2000 

3.5 

Helium/Argon 

2500 

2.9 

Helium/Argon 

3000 

2.5 

Hydrogen/Argon 
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Figure  6.  Terminal  velocity  of  A1?Q^  particles  in  argon 
as  a  function  of  diameter. 


2.3  Particle  Delivery  System 

The  dispersal  of  the  particles  within  the  shock  tube  test  section  is  an 
important  part  of  the  experimental  design  for  this  work.  It  is  by  no  means  a 
new  problem  and  previous  investigators  have  arrived  at  several  different  hard¬ 
ware  solutions  for  the  task.  These  solutions  include:  blowing  the  particles 
off  a  suspending  plate  in  a  shock  tube,  Calia  et  al.  (Ref.  8) ;  a  variety  of 
stirred  or  fluidized  bed  systems,  Cann  et  al.  (Ref.  11),  and  Lanzo  (Ref.  12); 
an  ultrasonic  levitating  system,  Harris  et  al.  (Ref.  13);  an  electrostatic 
levitating  system,  Adamo  et  al.  (Ref.  14);  air  jet  delivery  systems,  Decker 
et  al.  (Ref.  15);  and  Lowenstein  and  von  Rosenberg  (Ref.  16).  Each  of  these 
solutions  has  its  own  set  of  operating  peculiarities;  the  delivery  system  for 
this  work  is  based  upon  the  operating  principles  of  these  previously  designed 
systems  and  the  specialized  needs  of  the  shock  tube  experiment.  Specific  con¬ 
siderations  which  guided  the  particle  delivery  system  design  are:  1)  produc¬ 
tion  of  a  uniform  distribution  of  particles  in  the  viewing  area  for  optical 
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diagnostics,  2)  separation  of  shock  tube  and  particle  delivery  system  opera¬ 
tion  to  allow  the  delivery  system  development  to  proceed  independently  of  the 
shock  tube,  3)  a  system  design  which  was  fully  cognizant  of  the  possibility  oi 
particle  aggloraera tion  and  incorporated  previous  workers'  empirical  experience 
to  prevent  as  much  agglomera tion  as  possible,  4)  the  system  should  be  easy  to 
use  from  a  particle  handling  standpoint,  and,  finally,  5)  the  total  mass  of 
particles  required  for  one  experiment  should  be  substantially  less  than 
one  gram.  These  considerations  arise  from  a  detailed  examination  of  the  envi¬ 
ronmental  portion  of  the  experimental  design.  The  specific  solutions  to  these 
considerations  will  be  discussed  in  conjunction  with  the  final  particle 
delivery  system  design. 

A  radial  particle  injection  system  is  the  heart  of  the  particle  delivery 
system.  This  concept  is  unique  and  properly  capitalizes  on  shock  tube  geome¬ 
try  to  produce  a  uniform  particle  distribution  over  a  small  portion  of  the 
shock  tube  test  section  length.  In  contrast,  axial  delivery  systems  will 
produce  a  cloud  of  particles  that  is  unconstrained  in  length  and  diameter. 

For  example,  a  particle  jet  located  in  the  end  plate  of  the  test  section  will 
produce  a  relatively  long  cloud  of  particles.  This  cloud  will  follow  the 
spreading  of  the  jet,  which  forms  an  approximately  10°  half-angle  cone  accord¬ 
ing  to  Pai  (Ref.  17).  Figure  7  illustrates  this  type  of  delivery  system.  The 
path  length  through  the  particulate  cloud  varies  with  the  diametrical  path 
chosen  and  becomes  smaller  close  to  the  end  plate.  In  addition,  the  uncon¬ 
strained  length  of  the  particulate  cloud  would  dictate  larger  sample  loadings 
for  the  delivery  system.  These  considerations  point  to  a  radial  injector 
which  delivers  particles  over  a  small  axial  length  and  produces  a  uniform  dis¬ 
tribution  which  totally  fills  the  circular  cross  section  of  the  tube. 

The  radial  injection  concept  relies  on  an  appropriate  jet  duration  to 
produce  an  even  particle  distribution.  Jet  duration  is  conveniently  described 
in  terms  of  penetration  depth,  as  done  by  Cetagen  et  al.  (Ref.  18)  in  his 
plasma  jet  work.  The  concept  is  depicted  in  Figure  8.  The  penetration  depth, 
L* ,  is  defined  as  the  ratio  of  the  volume  of  gas  delivered  to  the  test  section 
to  the  orifice  area  of  the  jet.  Under-penetration,  as  shown  in  Figure  8(a), 
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PARTICLE 

CLOUD 


produces  an  inadequate  amou.it  of  mixing  in  the  test  volume.  Over  penetration , 
shown  in  Figure  8(b),  produces  a  steady  jet  with  sustained  vortices.  Although 
the  mixing  produced  by  this  system  is  very  desirable,  the  sustained  vortices 
could  produce  an  unacceptable  amount  of  size  sorting.  The  solution,  shown  in 
Figure  8(c),  is  a  penetration  depth  greater  than  tube  diameter  but  not  great 
enough  to  create  a  steady  jet  flow.  This  optimal  penetration  depth  will  pro¬ 
mote  mixing  due  to  the  jet  interaction  with  the  wall.  Optimal  penetration 
will  be  determined  empirically  and  will  be  on  the  order  in  a  few  tube 
diame  ters . 

Two  relevant  criteria  for  the  proper  radial  injector  performance  are  that 
the  particles  remain  entrained  in  the  gas  flow  and  that  a  minimal  amount  of 
size  sorting  occurs.  The  entrainment  criterion  is  satisfied  as  long  as  the 
local  gas  stream  velocity  is  greater  than  the  particle  terminal  velocity. 

Size  sorting  can  be  quantified  by  considering  particle  motion  in  a  circular 
path  similar  to  that  experienced  when  the  flow  changes  direction  at  the  tube 
wall.  The  forces  acting  on  the  particle  are  the  Stokes  drag,  due  to  relative 
motion  of  the  fluid,  and  the  force  from  acceleration  due  to  the  circular 
motion.  This  balance  of  forces  dictates  the  relationship  between  radial  and 
tangential  velocity  given  by  Eq.  (20). 


where : 

vr  -  radial  velocity 
Pp  -  particle  density 
a  -  particle  radius 
vt  -  tangential  velocity 
r*  -  turning  radius  of  the  flow 
M  -  gas  viscosity 
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Given  this  relationship  between  velocities,  size  sorting  can  be  quanti¬ 
fied.  The  radial  displacement  of  a  particle,  as  depicted  in  Figure  9  along 
with  a  simplified  vortex  for  the  flow,  over  the  flow  length  n  r*  is  given  by 
Eq.  (21). 


A 


where 


(21) 


A  -  normalized  radial  particle  movement 

6  -  radial  particle  movement 

Figure  10  illustrates  the  operating  region  for  the  radial  injector. 

Local  gas  velocities  should  always  be  greater  than  the  particle  terminal 
velocity,  i.e.,  to  the  right  of  the  terminal  velocity  curve.  To  avoid  size 
sorting  in  the  first  vortex,  the  gas  velocity  of  the  vortex  should  be  less 
than  that  for  A  =  10  percent,  which  leaves  a  significant  acceptable  onerating 
range  for  the  local  gas  velocity  of  the  first  vortex.  This  size  sorting 
analysis  assumed  a  steady,  well  established  flow.  The  jet  injection  system  is 
actually  a  transient  flow  and  the  flow  reversal  vortex  will  not  be  well  estab¬ 
lished,  which  implies  that  the  size  sorting  analysis  will  overpredict  size 
sorting  for  the  system.  In  addition,  the  most  severe  case  of  size  sorting  has 
been  examined;  the  jet  is  three-dimensional  in  nature  and  the  vortex  examined 
is  the  most  severe  one  that  will  be  encountered.  Therefore,  suitable  operat¬ 
ing  velocities  for  the  vortex  may  be  as  large  as  those  given  by  the 
A  =  10  percent  curve.  For  a  particle  10  pm  in  radius  this  implies  a  velocity 
in  the  vortex  of  approximately  30  cm/s.  Given  a  jet  orifice  size  and  Eq.  (22) 
from  Pai  (Ref.  17)  which  relates  the  centerline  velocity  of  the  jet  at  a  given 
axial  position  to  the  exit  velocity,  the  operating  parameters  for  the  jet  may 
be  determined.  They  are  listed  in  Table  6. 
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TEST  SECTION 


FIRST 

VORTEX 


Figure  9.  Relative  particle  motion  within  a  vortex. 


VELOCITY  (cm/s) 

Figure  10.  Operating  map  for  radial  injector,  calculation  of  A  assumes 
Pp  =  4  gm/cm3,  u  =  2.3  x  10~4  gm/cm~5,  and  r*  =  2  cm. 
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TABLE  6.  Operating  Parameters  for  the  Radial  Injection  System. 


Jet  Exit  Velocity 

-  82.5  cm/s 

Jet  Velocity  of  Vortex  Midpoint 

-  33  cm/s 

Jet  Orifice  Diameter 

0 . 6  cm 

Penetration  Depth 

-  30  cm  (3  tube  diameters) 

Delivery  Volume 

-  8.5  cm3 

Approximate  Delivery  Time 

-  0.36  s 

u 

c 

u 

e 


where : 


(22) 


uc  -  centerline  velocity  of  jet 
ue  -  velocity  of  jet  orifice 
d'  -  orifice  diameter 
x  -  axial  distance  from  the  orifice 


The  integration  of  the  radial  injector  to  the  shock  tube  experiment  is  a 
matter  of  timing;  the  experimental  events  must  be  sequenced  so  that  the  inci¬ 
dent  shock  arrives  at  the  optical  test  section  after  completion  of  particle 
delivery.  This  sequencing  will  not  be  difficult  because  of  the  relative  sta¬ 
bility  of  the  particle  dispersion  once  it  is  established.  The  stability  is 
directly  attributable  to  the  low  particle  terminal  velocities;  an  example  is  a 
10  pm  diameter  particle  of  AI2O3.  It  will  settle  one  centimeter  in  1.0s  when 
suspended  in  Argon  at  atmospheric  pressure  and  room  temperature.  Therefore, 
the  shock  would  have  to  arrive  at  the  particle  dispersion  within  200  ms  of  the 
dispersal  becoming  stable  and  stagnant. 
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The  proposed  particle  delivery  system  may  be  operated  independently  of 
the  shock  tube,  thus  reducing  the  set  up  time  for  each  injector  firing.  In 
addition  to  this  reducion  of  time  between  injector  firings  during  development, 
the  added  benefit  of  complete  optical  access  is  accomplished  by  firing  the 
injector  into  a  plexiglass  mock-up  of  the  shock  tube.  The  evolution  of  the 
particle  dispersal  can  therefore  be  monitored  and  the  results  used  to  fine- 
tune  the  operation  of  the  injector. 

The  jiarticle  delivery  system  has  been  defined  in  terms  of  the  character¬ 
istics  of  a  radial  jet  and  the  use  of  this  jet  within  a  shock  tube.  The  hard¬ 
ware  which  produces  the  jet  can  now  be  defined.  A  conceptual  drawing  of  the 
radial  injector  is  shown  in  Figure  11.  The  design  is  very  simple  in  concept. 
The  piston  is  driven  into  the  delivery  volume  producing  a  gas  jet  at  the  ori¬ 
fice  exit.  Particles  are  entrained  in  an  annular  volume  centered  on  the 
orifice  tube  and  through  which  50  percent  of  the  flow  passes.  This  annular 


PARTICLE  HOLDING 
VOLUME  AND 
ENTRAINMENT  SECTION 


Figure  1 1 .  Cross  section  of  the  particle  injector. 
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TABLE  7.  Radial  Injector  Specifications. 


Displacement  Volume 

(25-mm  bore,  nominal  stroke  17  mm, 

adjustable  for  variable  displacement  volume) 

8.5  cm ^ 

Orifice  Outlet  Diameter 

6  mm 

Annular  Volume 

(for  particles  to  be  injected) 

maximum  0.5  cm 3 

Annulus  Dimensions: 

Internal  Diameter 

6  mm 

Outside  Diameters 

Full  Depth 

Zero  Depth 

17  mm 

19  mm 

Orifice  Feed  Jet 

7.25  mm  diameter  screened 
with  #200  mesh.  Flow 

Area  -  14  mm2 

Annulus  Feed  Jets 

80  equally  spaced  on  a 
17.5-mm  diameter  circle, 
0.5-mm  diameter  jets. 

Total  Flow  Area  -  15.7  mm2 

Total  Jet  Inlet  Flow  Area 

29.7  mm2 

Total  Jet  Outlet  Flow  Area 

( 6-mm  diameter  orifice) 

28.3  mm2 

Approximate  Flow  Volumes  Through  the  Annulus 

10 

volume  can  accommodate  a  wide  range  of  particle  loadings  and,  if  necessary, 
could  be  modified  for  exceptionally  high  mass  loadings.  The  flow  enters  the 
annular  volume  through  80  equally  spaced  0.5-mm  holes.  To  facilitate  circum¬ 
ferential  flow  spreading,  the  jets  will  impinge  upon  a  ramp  in  the  annulus. 
Nominal  sizes  for  the  injector  are  given  in  Table  7. 


The  piston  will  be  driven  with  an  off-the-shelf  pneumatic  actuator.  Flow 
and  pressure  to  this  linear  air  motor  will  control  injection  times  and  hold 
the  jet  piston  in  place  during  each  shock  event. 
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The  radial  injector  will  mount  on  the  top  of  the  shock  tube;  the  annulus 
will  be  loaded  by  partially  disassembling  the  injector.  Particle  volumes 
required  for  the  experiments  are  dependent  upon  the  optical  properties  of  the 
particles  and  the  volume  into  which  they  are  delivered.  Based  on  the  angular 
spread  of  the  jet,  the  delivery  volume  will  be  approximately  1000  cm3  (1.5 
tube  diameters).  The  optical  depth  of  the  particle  dispersion  should  be 
approximately  0.1.  This  value  ensures  that  extinction  and  emission  measure¬ 
ments  are  linearly  dependent  on  the  response  of  each  particle  in  the  field  of 
view.  Optical  depths  significantly  less  than  0.1  will  result  in  an  overall 
decreased  signal  level;  greater  optical  depths  produce  a  nonlinearity  in  the 
optical  response  due  to  multiple  scattering  and  coincidence  effects.  This 
optical  depth  requirement  determines  the  number  density  as  previously  speci¬ 
fied  by  Eq.  (15).  The  particle  volumes  required  range  from  8  x  10“4  cm3  for 
0.1  pm  radii  particles  to  0.02  crn3  for  10.0  pm  radii  particle  (0.0032  gm  to 
0.16  gm  assuming  4  gm/cm3)  .  This  range  can  clearly  be  accommodated  in  the 
annulus. 

The  formation  of  particle  agglomerates  must  be  discouraged  by  the  parti¬ 
cle  handling  procedure  and  the  delivery  characteristics  of  the  injector.  The 
approach  is  to  take  all  the  necessary  steps  to  produce  an  agglomerate  free 
distribution.  Previous  workers'  experience,  Longacre  (Ref.  19),  has  shown 
that  water  significantly  increases  the  tendency  to  agglomerate.  To  combat 
this  effect  all  samples  will  be  dried  before  they  are  loaded  into  the  injec¬ 
tor.  Prior  to  filling  the  test  section,  the  test  section  and  injector  will  be 
evacuated.  The  particles  will  therefore  be  exposed  to  a  vacuum,  and  the 
injector  heated  to  remove  any  remaining  moisture.  The  delivery  system  is  con¬ 
figured  to  provide  a  turbulent  environment,  due  to  the  orifice  screen  and  the 
overall  turbulent  nature  of  the  jet.  This  will  further  serve  to  reduce 
agglomeration.  Finally,  the  particle  dispersal  will  be  subjected  to  the  inci¬ 
dent  and  reflected  shock  which  will  also  break  up  agglomerates. 

In  summary,  the  particle  injector  unit  is  simple,  both  in  concept  and  in 
terras  of  hardware  requirements.  It  addresses  all  of  the  relevant  constraints 
for  this  experiment.  The  design  is  based  upon  the  experience  of  other  workers 
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who  have  prouaced  dispersals  of  small  particulates  and  the  design  is  an 
adaptation  of  their  experience  to  the  specialized  needs  of  this  experiment. 
The  injector  can  be  developed  independently  of  the  experimental  shock  tube 
sequence,  it  capitalizes  on  shock  tube  geometry  to  produce  a  uniform  particle 
distribution;  it  requires  low  mass  loadings;  and  it  permits  significant 
precautions  to  eliminate  agglomeration. 
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3.  OPTICAL  DATA  ANALYSIS 


Mie  scattering  theory  defines  the  radiation  interaction  with  particles  in 
terms  of  the  complex  index  of  refraction  of  the  material.  This  formalization 
of  the  scattering,  absorption,  and  emission  process  is  well  documented  by  both 
van  de  Hulst  (Kef.  20)  and  Kerker  (Ref.  21),  and  the  mathematics  which 
describe  this  interaction  are  contained  in  a  computer  program  written  by  Dave 
(Ref.  22).  The  Mie  scattering  computations  provide  the  basic  ingredients  for 
this  experiment.  A  detailed  test  of  the  Mie  scattering  code  was  performed  by 
comparing  results  of  Wyatt  and  Stull  (Ref.  23).  The  agreement  in  predicted 
cross  sections  was  excellent;  all  results  were  equal  to  the  first  three  sig¬ 
nificant  digits.  In  addition,  the  predicted  radiation  intensity  as  a  function 
of  angle  was  verified  graphically  with  the  results  in  the  above  reference.  We 
feel  that  this  code  verification  is  exceptionally  important  since  both  the 
optical  design  and  index  of  refraction  determination  are  dependent  upon  the 
code . 


Estimated  index  of  refraction  data  is  required  to  predict  the  optical 
response  of  the  particulate  dispersion.  The  complex  index  of  refraction  is  a 
function  of  wavelength  and  temperature.  Particle  size  does  not  influence  the 
refractive  index.  Although  the  available  data  is  sketchy  for  rocket  exhaust 
particulates,  a  reasonable  guess  is  provided  by  the  following  references: 

Adams  (Ref.  14),  Friese  (Ref.  24),  Dalzell  and  Sarofim  (Ref.  25),  Lee  and  Tien 
(Ref.  26),  Ludwig  et  al.  (Ref.  27),  Whitson  (Ref.  28),  Mularz  (Ref.  3),  and 
Pluchino  et  al.  (Ref.  2).  Figure  12  illustrates  the  index  of  refraction  data 
available  from  these  references  for  AI2O3  as  a  function  of  temperature  and 
wavelength.  In  several  cases  only  one  part  of  the  two  refractive  indices 
(n,k)  is  available,  and  data  from  different  sources  are  not  entirely  consis¬ 
tent.  A  similar  amount  of  information  is  available  for  carbon;  a  very  limited 
amount  of  data  is  available  for  MgO  and  ZrC>2,  and  no  data  has  been  found  for 
ZrC.  The  available  data  has  been  used  to  construct  several  data  sets  which 
were  then  used  to  predict  the  optical  response  of  particulate  dispersions.  A 
typical  set  of  index  of  refraction  curves  is  given  for  AI2O3  at  1000  K  in 
Figures  13  and  14. 
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Figure  12.  A1?0^  refractive  index  data  availability. 


REAL 

INDEX 


0.1  0.3  1.0  10.0 


WAVELENGTH  (pm) 


Figure  13.  Real  index  of  refraction  range  for  Al2°i  °f  1000  K 
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Figure  14.  Imaginary  index  of  refraction  range  for  A1?Q^  of  1000  K. 

Mie  scattering  theory  is  used  to  predict  the  optical  response  of  a  parti¬ 
cle  or  set  of  particles  given  the  index  of  refraction,  particle  size,  and 
radiation  wavelength.  Particle  size  and  radiation  wavelength  are  combined 
into  one  quantity,  the  size  parameter.  Given  these  two  inputs,  refractive 
index  and  size  parameter,  the  optical  response  of  a  single  particle  is  yiven. 
This  response  is  represented  by  the  extinction  efficiency,  ye ,  the  absorption 
efficiency,  Qg ,  and  the  angular  scattering  distribution  of  radiation  for  the 
two  polarization  states.  The  extinction  of  radiation  through  a  particulate 
dispersion  is  dependent  on  extinction  efficiency,  the  number  density  of  parti¬ 
cles,  and  the  optical  path  length,  as  given  by  the  Beer-Lambert  Law  formalism 
in  Eq.  (23).  This  equation  assumes  all  particles  to  be  the  same  size;  if  a 
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size  distribution  is  present,  the  exponential  operand  is  replaced  with  an 
integral  which  represents  the  size  distribution,  Eq.  (24). 

7—  =  exp( -Q  N  ira2£)  (23) 

I  e 

o 

where 

I  -  observed  radiation  intensity 
IQ  -  radiation  intensity  when  no  extinction  occurs 
Qe  -  extinction  efficiency 
N  -  number  density  of  particles 
a  -  particle  radius 
l  -  optical  path  length. 

—  «  exp  (  /  °°  -  QeU)  N'(a)  ia2f  da)  (24) 

o 

where  N' (a)  represents  the  size  distribution  of  particles:  N  »  f N' (a)  da. 
The  effects  of  a  size  distribution  will  be  discussed  later  in  this  section; 
until  then  a  monodisperse  system  will  be  assumed. 

The  absorption  efficiency,  is  directly  related  to  the  emission  from 

the  particle  dispersion.  Kirchoff's  law  states  that  the  spectral  absorptivity 
is  equal  to  the  spectral  emissivity.  Therefore,  the  spectral  emissivity  can 
be  determined  from  the  absorption  efficiency,  as  given  by  Eq.  (25). 

2 

e  =  1  -  exp(-Nna  Q  l )  (25) 

A  a 

where 


-  spectral  emissivity 
Qa  -  absorption  efficiency 

Scattering  is  also  predicted  by  Mie  theory.  This  prediction  takes  the 
form  of  an  angular  light  distribution  for  the  two  perpendicular  polarization 
components  of  light.  This  implies  that  relating  a  scattered  light  measurement 
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to  the  results  of  a  Mie  theory  calculation  is  potentially  more  difficult. 
Extinction  and  emission  measurements,  related  to  Mie  theory  calculations  by  Qg 
and  Qa  ,  are  independent  of  polarization  state  and  angular  position.  Scat¬ 
tering  measurements  are  dependent  on  polarization  and  angle;  therefore  a 
scattering-based  measurement  is  potentially  more  difficult  to  use  for  compar¬ 
ing  measurements  to  calculation  results. 

The  three  basic  measurement  types--extinction ,  emission,  and  scattering  — 
can  be  compared  with  the  results  of  Mie  theory  calculations  and  have  been 
introduced.  In  the  remaining  portion  of  this  section,  a  method  for  inferring 
the  complex  index  of  refraction  from  the  results  of  these  measurements  will  be 
developed. 


3.1  Sensitivity  of  Optical  Measurements  to  the 
Index  of  Refraction 

The  overall  goal  for  this  project  is  to  produce  an  experimental  program 
that  will  provide  as  direct  a  measure  as  possible  of  the  indices  of  refraction 
for  small  particles.  The  next  step  in  this  development  is  to  fully  understand 
the  sensitivity  of  the  optical  measurements  to  the  indices  of  refraction. 

3.1.1  Extinction  and  Emission  Measurements  -  A  range  of  emission  and  extinc¬ 
tion  responses  is  given  in  Figures  15  and  16.  These  curves  were  generated 
based  on  the  indices  of  refraction  shown  in  Figures  13  and  14.  The  absorption 
efficiency  reflects  the  large  range  of  possible  imaginary  indices  having  a 
large  range  of  possible  values.  The  extinction  efficiency,  which  includes 
scattering  effects,  has  a  much  smaller  response  range,  indicative  of  the 
smaller  range  in  the  real  index.  When  examining  the  sensitivity  of  optical 
response  to  the  index  of  refraction,  remember  that  there  are  two  variables 
that  determine  the  optical  response:  size  parameter  and  complex  index  of 
refraction.  The  size  parameter  can  influence  the  optical  response  to  a  very 
large  degree.  A  mapping  of  the  size  parameter  as  a  function  of  particle 
radius  and  wavelength,  shown  in  Figure  17,  graphically  illustrates  that  the 
range  of  possible  size  parameters  is  from  0.1  to  100.  As  an  input  to  the  Mie 
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Figure  15.  Absorptive  efficiency  for  Al?Oi  at  1000  K. 


Figure  16.  Extinction  efficiency  for  Al?Q->  at  1000  K. 
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X  (urn) 

Figure  17.  Hie  scattering  size  parameters,  a ,  as  a  function 
of  particle  radius  and  wavelength. 

calculations  the  size  parameter  combines  the  effect  of  particle  radius  and 
radiation  wavelength  into  one  parameter,  but  obscures  the  fact  that  the  index 
of  refraction  is  a  function  of  wavelength  only.  For  a  given  wavelength,  the 
optical  response  of  a  system  can  vary  due  to  changes  in  particle  size  (i.e., 
size  parameter)  only.  This  fact  suggests  that  an  optical  measurement's 
response  could  be  influenced  by  the  size  parameter  alone.  This  variation  of 
response  with  particle  radius  (size  parameter)  can  be  used  to  infer  specific 
information  about  the  refractive  index  at  a  given  wavelength. 


The  range  of  possible  extinction  efficiencies  is  given  in  Figure  18. 

Four  indices  of  refraction  representative  of  a  large  range  were  used  to  gener¬ 
ate  this  figure.  The  small  imaginary  indices  represent  of  nonabsorbing  mater¬ 
ials  and  the  large  imaginary  indices  represent  absorbing  materials.  The  two 
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Figure  18.  Extinction  efficiencies  for  a  wide  range  in  indices 
of  refraction. 


real  components  for  the  refractive  index  represent  the  range  typically 
observed  for  the  real  index.  The  response  range  is  exceptionally  large  at 
small  size  parameters  but  narrows  as  the  size  parameter  becomes  larger.  This 
apparent  insensitivity  of  the  extinction  to  the  refractive  index  at  large 
values  of  the  size  parameter  is  classically  known  as  the  geometric  limit.  The 
insensitivity  of  the  extinction  efficiency  to  the  optical  properties  allows 
extinction  data  for  large  size  parameters  to  be  used  as  a  direct  measurement 
of  the  product  of  number  density  and  particle  cross  sectional  area,  Eq.  (26). 

ln( I/I  ) 

Nira^i  =  { - - — (26) 

e 

The  sensitivity  of  the  extinction  measurement  to  refractive  index  is 
further  explored  in  Figures  19  -  23.  These  figures  represent  slices  of  the 
response  surface  at  constant  imaginary  or  real  index  for  several  size  param¬ 
eters.  Extinction  efficiency  is  demonstrated  to  be  a  function  of  both  the 
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Figure  23.  Extinction  efficiency  for  imaginary  index  of  2,0  and 
variable  real  index. 


real  and  imaginary  indices;  it  is  relatively  insensitive  to  either  index  for  a 
size  parameter  equal  to  100.  This  further  supports  the  use  of  large  size 
parameter  measurements  to  determine  the  number  density-area  product.  The 
largest  variation  in  response  for  a  variable  imaginary  index,  as  shown  by 
Figures  19  and  20,  occurs  at  a  size  parameter  equal  to  one.  A  decreased  sen¬ 
sitivity  for  size  parameters  other  than  one  is  apparent  for  larger  real  index 
values.  Figures  21  through  23  illustrate  the  response  range  for  variable  real 
index  at  several  size  parameters;  once  again,  the  response  is  insensitive  to 
the  index  of  refraction  at  large  size  parameters.  The  responses  shown  in 
these  plots  also  illustrate  a  decrease  in  dynamic  range  as  the  material 
becomes  more  absorbing  and  as  with  the  previous  set  of  plots,  the  highest  sen¬ 
sitivity  to  the  real  refractive  index  occurs  for  a  size  parameter  equal  to 
one.  These  sensitivity  plots  show  that  the  extinction  efficiency  can  be  used 
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as  a  measure  of  the  number  density  cross-sectional  area  product  for  large  size 
parameters  and  is  a  measure  of  the  real  refractive  index  for  size  parameters 
equal  to  one. 

A  special  note  must  be  made  regarding  the  extinction  response  for  a  non¬ 
absorbing  material  as  represented  in  Figure  21.  The  response  range  for  small 
values  of  the  real  index  and  size  parameter  equal  to  100  is  larger  than  demon¬ 
strated  In  any  of  the  other  cases  and  could  be  a  error  source  in  a  measurement 
of  the  number  density-area  product.  Fortunately  this  combination  of  refrac¬ 
tive  indices  is  physically  unlikely  for  a  large  size  parameter  with  rocket 
exhaust  particles.  A  size  parameter  of  100  only  occurs  for  near  ultraviolet 
wavelengths  with  our  particle  size  range;  the  real  refractive  index  in  this 
range  is  normally  greater  than  1.4.  The  large  variation  in  response  occurs 
between  1.05  and  1.2;  therefore,  extinction  is  still  a  reasonable  measure  of 
the  number  density-area  product  for  a  nonabsorbing  material.  In  an  actual 
experimental  system,  the  final  indices  that  are  determined  would  be  used  to 
verify  this  assumption. 

A  similar  study  of  absorption  efficiencies  and  their  dependence  on  size 
parameter  and  refractive  index  is  given  in  Figures  24  through  27.  The  case  of 
a  zero  imaginary  refractive  index  has  been  omitted  since  a  zero  imaginary 
index  indicates  a  material  with  no  absorptive  properties.  The  first  two  plots 
show  the  variation  of  absorption  efficiency  with  imaginary  index.  Once  again 
the  most  dramatic  response  range  is  found  for  a  size  parameter  equal  to  one. 
However,  there  is  no  area  of  insensitivity,  as  there  was  with  extinction. 

This  indicates  that  emission  measurements  over  the  entire  range  of  refractive 
index  will  yield  information  about  the  imaginary  index.  Figures  26  and  27 
show  the  response  of  the  absorption  efficiency  to  the  real  index.  It  is  not 
surprising  to  find  that  the  emission  response  range,  represented  by  the 
absorptive  efficiency,  is  not  as  dramatic  as  for  the  imaginary  index.  In  sum¬ 
mary,  emission  measurements  will  be  primarily  dependent  on  the  imaginary  index 
and  secondarily  dependent  on  the  real  index. 
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Figure  26.  Absorption  efficiency  for  imaginary  index  of  0.50 
and  variable  real  index. 


Figure  27.  Absorption  efficiency  for  imaginary  index  2.0  and 
variable  real  index. 
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The  previous  diaeiisuion  luin  ox.iinl nod  tin?  uunulLlvlLy  ol  the  extinct  ion 
and  emissive  response  of  systems  of  particles  to  the  refractive  index  of  the 
materials.  The  conclusions  from  this  study  are  summarized  in  Figure  28  which 
presents  extinction  and  emission  sensitivities  grouped  by  size  parameter 
range.  For  large  size  parameter  values,  the  extinction  measurement  can  be 
used  to  determine  the  number  density-area  product  for  the  system  and  the  emis¬ 
sion  measurement  is  a  reasonable  measure  of  the  imaginary  refractive  index. 
When  10  >  a  >  1,  the  extinction  measurement  becomes  sensitive  to  the  real  and 
imaginary  indices,  and  the  emission  measurement  is  sensitive  to  the  imaginary 
and  real  indices.  When  a  <  1,  both  measurements  are  sensitive  to  the 
imaginary  and  real  indices. 
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Figure  28.  Extinction  and  emission  measurements  and 
their  sensitivities. 
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3.1.2  Scattering  Measurements  -  The  scattering  response  is  a  function  of 
angle  and  polarization  in  addition  to  size  parameter  and  index  of  refraction. 
Due  to  the  complexity  of  the  scattering  response,  several  measurement  types 
must  be  considered  along  with  their  sensitivities  to  the  refractive  index. 
These  measurements  may  be  compactly  represented  by  the  polarization  ratio  and 
the  angular  intensity  of  scattered  natural  light. 


Figures  29  through  32  display  the  polarization  ratio  as  a  function  of 
angle  for  four  refractive  indices  at  each  of  four  size  parameters.  Polariza¬ 
tion  ratio  is  defined  in  Eq.  (27). 
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where 


P  -  polarization  ratio 
1^  -  perpendicular  intensity 
I||  -  parallel  intensity 

Figures  29  and  30  show  an  insensitivity  of  response  to  the  refractive  index 
for  size  0.1  and  1.0.  At  larger  size  parameters  the  response  becomes  depen¬ 
dent  on  the  refractive  index  in  addition  to  becoming  a  highly  structured  func¬ 
tion  of  angle,  as  shown  in  Figures  31  and  32.  The  degree  of  structure  in 
these  responses  implies  that  although  in  principle  refractive  index  influences 
the  response,  measurements  based  on  these  predictions  would  be  very  prone  to 
error  due  to  the  large  response  changes  for  small  changes  in  angle.  The 
structure  becomes  more  detailed  for  a  size  parameter  of  100,  but  the  response 
becomes  sensitive  to  the  imaginary  index. 

The  shift  in  the  function  from  a  size  parameter  of  1.0  to  10.0  also  sug¬ 
gests  that  the  measurement  could  be  used  to  diagnose  the  size  of  the  particles 
in  the  dispersion.  However,  this  measurement  would  be  useful  only  in  the  size 
range  where  the  response  function  is  making  the  transition  from  the  invariant 
curve  for  a  size  parameter  of  one  to  the  robust  function  for  size  parameter 
equal  to  ten. 
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The  possible  uses  for  scattering  polarization  measurements  are  summarized 
in  Figure  33.  This  shows  that  the  measurement  is  insensitive  to  particle  size 
or  refractive  index  for  size  parameters  less  than  1.0.  For  size  parameter 
between  1.0  and  10.0,  the  measurement  could  be  used  for  sizing  or  as  a  measure 
of  refractive  index,  but  it  is  exceptionally  sensitive  to  angular  position. 

At  size  parameters  between  10  and  100  the  responses  become  sensitive  to  the 
imaginary  index  but  are  still  dependent  on  angular  position. 
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Figure  33.  Polarization  sensitivity. 
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Scattered  intensity  of  incident  natural  light  forms  the  other  type  of 
scattering  measurement  that  can  be  made.  These  measurements  can  be  either 
absolute  intensity  determinations  at  a  specific  angle  or  a  ratio  of  intensi¬ 
ties  at  two  or  more  angles.  The  latter  measurement  type  has  the  advantage  of 
being  independent  of  the  number  density-area  product  for  the  particle  disper¬ 
sion.  Plots  of  intensity  as  a  function  of  angle  for  four  indices  of  refrac¬ 
tion  are  given  in  Figures  34  through  37  for  size  parameters  0.1,  1.0,  10.0  and 

100.  The  intensity  distribution  is  a  well  behaved  function  for  size  param¬ 
eters  0.1  and  1.0.  In  addition,  the  influence  of  the  two  refractive  indices 
on  the  response  is  apparent.  The  real  index  dramatically  affects  the  magni¬ 
tude  of  the  response  and  the  imaginary  index  has  almost  no  effect  on  the 
response  at  all.  The  curve  shape,  commonly  called  the  dissymmetry  ratio,  is 
represented  by  the  ratio  of  intensities  at  two  angles  and  is  insensitive  to 
the  refractive  index.  However,  this  ratio  is  reasonably  sensitive  to  the  size 
parameter  and  could  therefore  be  used  as  a  diagnostic  tool  for  particle  size 
at  a  fixed  wavelength.  For  a  size  parameters  of  10.0,  represented  by 
Figure  36,  the  uses  for  scattering  measurements  become  much  less  clear;  the 
angular  intensity  function  becomes  much  more  structured  and  its  use  as  a  diag¬ 
nostic  aid  becomes  more  tenuous.  As  with  the  polarization  ratio,  the  highly 
structured  curve,  although  sensitive  to  the  refractive  index  in  some  regions, 
is  difficult  to  use  for  determinative  purposes.  As  an  example,  the  first 
scattering  lobe  for  the  size  parameter  equal  to  10.0  is  reasonably  sensitive 
to  the  real  refractive  index.  However,  any  intensity  determination  at  a 
specific  angle  would  require  a  very  high  degree  of  accuracy  in  specifying  the 
exact  angle  for  the  scattering  because  of  the  function's  large  slope  in  this 
region.  Therefore,  the  scattering  measurement  is  not  a  good  diagnostic  aid 
for  the  refractive  index  at  high  size  parameter  values. 

Finally,  the  intensity  distribution  for  a  size  parameter  of  100.0  is 
shown  in  Figure  37.  The  salient  features  of  this  plot  are  a  fine  scale  struc¬ 
ture  superimposed  on  a  larger  scale  function  variation.  The  response  sensi¬ 
tivity  for  this  size  parameter  has  shifted  to  the  imaginary  index  in  the  first 
ninety  degrees  of  the  function. 
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Scattering  measurements  and  their  utility  are  summarized  in  Figure  30. 

For  size  parameters  between  0.1  and  1.0  the  intensity  function  is  well  behaved 
and  sensitive  to  the  real  index.  The  shape  of  the  curve  is  sensitive  to  the 
size  parameter.  For  size  parameters  between  1.0  and  10.0,  a  transition  occurs 
to  a  structured  function  that  is  sensitive  to  the  real  index  in  the  first  Lolx' 
of  the  function.  At  larger  size  parameters,  the  scattering  becomes  even  more 
structured  and  the  sensitivity  shifts  to  the  imaginary  index  in  the  first 
ninety  degrees  of  scattering. 
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Figure  38.  Scattering  intensity  sensitivities. 
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3.2  Data  Inversion:  Determination  of  the  Refractive  Index 


Tl'.  measurements  discussed  in  the  previous  section  can  be  used  to  deter¬ 
mine  the  refractive  index  of  the  particles  in  a  dispersion.  In  principle, 
three  unknowns  will  be  present  for  any  system  of  particles,  the  two  components 
of  the  refractive  index  and  the  number  density-area  product.  Since  the 
refractive  index  is  dependent  on  the  wavelength  only  and  not  the  particle 
size,  several  experiments  with  the  same  material  but  different  size  particles 
may  be  used  to  for  data  inversion.  The  results  of  several  experiments  are 
thus  used  in  concert  to  determine  the  refractive  index. 

The  most  direct  measurements,  as  presented  in  the  previous  section,  are 
extinction  and  emission.  In  general,  these  two  measurements  will  form  the 
backbone  of  the  experiment.  Where  necessary,  scattering  measurements  can  be 
used  to  fill  data  gaps.  The  sensitivity  variation  with  size  parameter  can  be 
used  to  good  advantage.  The  key  to  determining  of  the  refractive  index  is  to 
use  the  extinction  as  a  measure  of  number  density-area  product  where  it  is 
insensitive  to  the  refractive  index.  This  will  be  in  the  near  ultraviolet  for 
large  particles  as  shown  in  Figure  39.  Given  this  absolute  determination  of 
the  number  density-area  product,  emission  and  extinction  measurements  at  other 
wavelengths  are  indicative  of  the  refractive  indices  at  these  wavelengths.  In 
a  sense  the  experiment  is  based  on  the  transition  of  sensitivity  with  decreas¬ 
ing  size  parameter.  The  results  of  experiments  at  smaller  particle  sizes  are 
then  used  in  conjunction  with  the  large  particle  experiment  to  determine  the 
refractive  index  over  a  large  wavelength  range. 

An  indicative  example  is  given  by  a  two  experiment  system:  one  for  par¬ 
ticle  sizes  of  10  pm  and  one  for  particle  sizes  of  0.3  pm.  For  the  10  pm 
particle  size  experiment,  the  short  wavelength  extinction  measurements  give 
the  number  density  area  product.  At  longer  wavelengths  the  extinction  and 
emission  measurements  can  be  used  to  infer  the  refractive  index.  This  proced¬ 
ure  of  inferring  the  indices  of  refraction  based  on  the  optical  response  uses 
the  Mie  calculations  and  is  subject  to  the  sensitivities  associated  with  the 
size  parameter;  as  the  size  parameter  decreases,  the  results  are  wore 
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Figure  39.  Multiple  experiment  determination  of  refractive  index. 

sensitive  to  the  refractive  index.  At  10  ym  the  emission  measurement  is  rea¬ 
sonably  sensitive  to  the  imaginary  index  and  the  extinction  measurement  yields 
a  range  for  the  real  index.  These  results  are  used  in  conjunction  with  those 
of  the  0.3  ym  particle  experiment.  The  refractive  index  range  for  a  10  ym 
wavelength,  from  the  10  ym  particle  experiment  and  the  extinction  and  emission 
measurement  at  10  ym  from  the  0.3  ym  particle  experiment  yields  a  range  for 
the  number  density-area  product  for  the  dispersion  of  0.3  ym  particles.  Thus, 
the  results  from  the  10  ym  particle  experiments  are  used  to  tie  down  the  num¬ 
ber  density-area  product  for  the  0.3  ym  dispersion.  Measurements  at  other 
wavelengths  are  used  to  infer  the  refractive  index  at  those  wavelengths  and 
the  results  of  the  two  experiments  intersect  to  produce  a  smaller  range  of 
possible  indices.  Additional  experiments  at  other  particle  sizes  can  then  be 
done.  One  of  two  methods  is  used  to  determine  the  number  density-area  product 
for  the  system:  either  the  extinction  at  short  wavelengths  or  the  refractive 
index  at  long  wavelengths.  The  method  which  provides  the  greater  accuracy  is 
the  one  which  is  chosen. 
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As  a  fur  the  r  illustration  o£  the  multiple  experiment,  da  La  inversion  proc¬ 
ess,  two  simulated  inversions  have  been  performed:  one  tor  alumina  and  one 
for  carbon.  The  index  of  refraction  data  base  which  was  previously  discussed 
was  used  to  produce  mock  experimental  data.  Appropriate  experimental  error 
bars  were  added  to  these  numbers,  and  an  inversion  was  performed  within  the 
limits  of  experimental  uncertainty  (error  bars  were  ±1  percent  for  extinction 
measurements  and  ±7  percent  for  emission  measurements,  these  errors  are  indic¬ 
ative  of  a  very  carefully  performed  experiment).  The  results  are  shown  in 
Figure  40  for  four  experiments  for  each  system,  10.0  pm,  3.0  pm,  1.0  pm,  and 
0.7  pm.  The  solution  for  each  system  is  in  the  form  of  a  solution  region. 

The  refractive  indices  within  these  regions  represent  an  extinction  and  emis¬ 
sion  response  range  that  is  of  the  same  order  as  the  extinction  and  emission 
errors  for  the  most  sensitive  experiment.  Each  of  these  results,  for  the 
refractive  index  at  10  pm  is  representative  of  an  extinction  and  emission 
response  range  of  approximately  ±1  percent  and  +7  percent  plus  the  error 
involved  in  the  number  density  area  product,  typically  1-5  percent.  It  is 
interesting  to  note  the  decrease  in  sensitivity  for  high  values  of  the  real 
index.  This  is  exemplified  by  the  much  larger  range  in  imaginary  and  real 
index  for  carbon  and  was  previously  mentioned  as  a  part  of  the  sensitivity 
study  portion  of  this  report.  However,  each  of  these  solution  regions  is 
representative  of  the  same  range  in  extinction  and  emission  response  of  a 
system  of  particles. 

Finally,  it  must  be  pointed  out  that  the  process  of  inverting  the  data  is 
a  two-dimensional  iteration  process.  The  extinction  and  emission  are  depen¬ 
dent  on  both  refractive  indices,  and  the  solution  process,  from  an  extinction 
or  emission  response  to  the  refractive  index,  is  iterative  in  nature.  For  a 
large  experimental  program  this  data  reduction  process  would  have  to  be  auto¬ 
mated;  a  two-dimensional  bisection  solution  process  as  used  by  Parker 
(Ref.  29)  would  be  appropriate. 
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Figure  40.  Data  inversion  results  for  the  refractive  index. 


3.3  Particle  Size  and  Shape  Effects 

The  analysis  so  far  in  this  report  has  ignored  the  effects  of  disperse 
particle  size  distributions  and  particle  shape.  The  effect  of  a  particle  size 
distribution  is  easily  accounted  for  while  the  effect  due  to  particle  shape  is 
a  little  more  tenuous.  In  general,  a  particle  size  distribution  merely  repre¬ 
sents  a  complication  in  the  calculation  procedure.  The  optical  response  must 
now  be  considered  to  be  that  of  a  sum  of  particles  which  have  a  range  in  radii 
and  efficiency  corresponding  to  the  size  distribution  of  the  particles.  If 
the  system  is  optically  thin,  the  response  is  simply  the  sum  of  the  responses 
for  the  varying  size  particles,  i.e.,  an  integral  of  emission  or  extinction 
over  the  size  distribution  is  appropriate.  The  important  aspect  of  size  dis¬ 
tribution  effects  is  knowing  the  appropriate  size  distribution  for  the  system. 
An  unknown  size  distribution  can  catastrophically  affect  the  error  in  the 
experiment. 

Figures  41  through  45  illustrate  the  effect  of  a  size  distribution  on 
the  extinction  and  absorbtion  efficiencies  for  a  particle  system.  Figure  41 
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4 1 .  Nominal  1  um  size  distribution  and  agglomerated  distribution 


Figure  42.  Extinction  efficiency  for  a  monodisperse  and 
size  distribution  type  systems. 


Figure  45.  Agglomeration  effect  on  absorption. 

shows  a  nominal  1  ym  distribution  (as  opposed  to  a  monodisperse  1  ym  system). 
Figures  42  and  43  show  the  effect  of  the  size  distribution  on  the  area 
weighted  average  extinction  and  absorption  efficiency;  the  extinction  effi¬ 
ciency  is  exceptionally  sensitive  to  the  change  from  a  monodisperse  system  to 
a  size  distribution.  An  agglomerated  system  has  also  been  modeled.  This  sys¬ 
tem  was  produced  by  assuming  that  the  total  number  of  particles  was  reduced  by 
a  factor  of  2.  The  size  was  accordingly  shifted  upward  by  the  cube  root  of  2. 
This  is  not  an  extreme  agglomeration  model.  In  som->  systems  there  is  a  ten¬ 
dency  to  agglomerate  to  one  specific  size.  The  effects  of  the  agglomeration 
are  shown  in  Figures  44  and  45.  Once  again  the  extinction  efficiency  is  rea¬ 
sonably  sensitive  to  the  change  in  size  distribution  while  absorption  effi¬ 
ciency  is  almost  unaffected.  An  additional  error  from  agglomeration  is  a 
possible  change  in  shape. 
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The  solution  to  the  error  caused  by  particle  size  distributions  is  to 
always  know  the  distribution  of  particles  in  the  experimental  system.  This 
knowledge  can  be  obtained  with  pre-  and  post-experiment  SEM  particle  studies, 
and  by  using  scattering  dissymmetry  as  a  diagnostic  indicative  of  size.  The 
size  distribution  of  the  particulates  can  be  improved  toward  the  ideal  of  a 
monodispersion  by  pre-sorting  the  particles.  This  sorting  could  be  performed 
in  a  cascade  irapactor  due  to  the  small  sample  quantities  required. 

Particle  shape  effects  are  much  harder  to  quantify;  strict  Mie  theory  as¬ 
sumes  perfect  spherical  particles  and  the  optical  response  deviation  due  to 
particle  shape  is  not  well  quantified.  However,  recent  work  by  Pinnick  et  al. 
(Ref.  30),  has  experimentally  determined  for  pure  scatterers  with  a  real  index 
of  1.54  that  random  particle  shape  has  very  little  effect  on  the  optical 
response  for  systems  of  particles  except  for  scattering  at  angles  between  30 
and  180  degrees  with  size  parameters  greater  than  5.  Therefore,  our  extinc¬ 
tion  and  emission  measurements  will  be  largely  unaffected  by  particle  shape. 


3.4  Measurement  Recommendations 


The  sensitivity  analysis  and  data  inversion  study  have  identified  which 
measurements  can  be  used  in  which  size  regimes  to  produce  information  relevant 
to  the  refractive  index.  This  information  is  summarized  in  Figure  40.  Note 
that  in  all  size  regimes  the  two  primary  measurements  are  extinction  and  emis¬ 
sion;  the  quantity  to  which  the  measurement  is  sensitive  may  change,  but  for 
all  sizes  these  two  measurements  are  useful.  Scattering  can  be  used  in  the 
small  size  parameter  range  as  a  diagnostic  aid  for  real  refractive  index  and 
size.  It  is  theoretically  useful  in  the  larger  size  regimes,  but  the  degree 
of  structure  in  the  angular  distributions  would  make  data  analysis  inversion 
and  analysis  extremely  error  prone.  Given  these  trends,  the  measurements 
listed  in  Table  8  are  recommended,  in  order  of  preference,  as  the  primary 
diagnostics  for  this  refractive  index  study. 
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Figure  46.  Optical  measurement  sensitivity 


TABLE  8.  Measurement  recommendations. 


• 

Extinction  => 

Ntt  a  2  £. 

n,  k 

• 

Emission  => 

n,k 

• 

Scattering  Assymetry  => 

size  information  for  a  <  1.0 

• 

Absolute  Scattering  > 

real  index  0.6  <  u  <  6.0 
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4.  DIAGNOSTIC  MEASUREMENT  TECHNIQUES 


4.1  Requirements  and  Proposed  Approach 

From  the  preceding  discussion  it  is  clear  that  diagnostic  techniques  must 
be  developed  in  two  areas.  The  primary  requirement  is  the  capability  to 
simultaneously  acquire  extinction,  emission,  and  scattering  measurements  over 
a  broad  wavelength  range,  i.e.,  0.1  pm  <  X  <  10  pm,  within  a  time  duration 
determined  by  a  shock  tube  test  time.  Secondly,  techniques  are  needed  for 
characterizing  aerosol  number  density  and  particle  size.  While  information  on 
number  density  and  particle  size  can  be  unfolded  from  the  extinction  and  scat¬ 
tering  data,  it  is  useful  to  consider  additional  techniques  that  can  provide 
such  information  independent  of  a  particle  optical  properties. 

The  approach  we  have  chosen  to  meet  these  objectives  is  to  design  a  shock 
tube  experiment  implemented  with  numerous  spectroscopic  and  radiometric  meas¬ 
urements  that  can  be  deployed  simultaneously  on  each  shot.  A  discussion  of 
the  proposed  measurement  techniques  is  presented  below. 

4.2  Description  of  Proposed  Diagnostic  Techniques 

4.2.1  Particle  Number  Density  and  Size  Characterization  -  In  the  previous 
chapter,  it  was  shown  that  with  the  appropriate  choice  of  particle  size  param¬ 
eter,  a,  (where  a  =  2ira/A  )  extinction  and  selected  scattering  measurements  can 
be  used  to  infer  information  on  particle  number  density  and,  to  a  more  limited 
extent,  check  assumptions  on  particle  size  distribution.  For  example,  at 
large  size  parameters  (u  >  10)  extinction  efficiencies  become  independent  of 
particle  size  and,  to  a  great  extent,  of  refractive  index.  Thus,  measured 
extinction  signals  can  be  useful  for  determining  line-of-sight  number  densi¬ 
ties.  At  intermedia te-to-small  size  parameters,  angular  scattering  distribu¬ 
tions  can  provide  information  on  particle  size  (see  Figures  34  through  37). 
Thus,  to  the  extent  it  is  practical,  particle  size  and  wavelength  should  be 
varied  in  the  experiments  and  analysis  of  the  observed  extinction  and  scatter¬ 
ing  signals  could  be  used  to  infer  particle  number  density  and  size 
information. 
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A  second  technique  which  should  be  used  to  characterize  particle  size 
distribution  is  to  perform  micrographic  particle  sampling/sizing  both  prior  to 
particle  loading  and  subsequent  to  shock  heating.  Scanning  electron  micro¬ 
graphy  in  addition  to  optical  micrography  is  recommended.  Post-test  analysis 
of  the  shock  heated  particulates  is  useful  for  establishing  limits  on  possible 
effects  due  to  particle  agglomeration  and/or  particle  fragmentation  in  the 
shock  tube. 

Finally,  it  is  recommended  that  high  resolution  snapshots  of  the  particle 
field  behind  the  shock  also  be  taken.  An  attractive  technique  is  to  take 
short  exposure  (TeXpOSure  <  20  ns)  holograms  using  a  commercially  available 
laser-based  holography  system.  The  short  exposure  time  made  possible  by  the 
laser  overcomes  the  problem  of  blurring  due  to  particle  motion.  The  laser 
also  provides  a  large  coherence  length  which  helps  to  maximize  the  depth  of 
field  for  the  holograms.  Thus,  using  this  technique  with  appropriate  recon¬ 
struction  of  the  holograms,  particle  size  and  number  density  information  can 
be  acquired  as  a  function  of  shock  tube  radius,  i.e.,  at  selected  planes  along 
the  line  of  sight.  Of  course,  this  three-dimensional  imaging  technique  can 
also  yield  information  on  particle  shape.  The  holographic  technique  is  most 
useful  for  particle  radii  >  3  u™  and  scattering  optical  depths  <  0.1  (negli¬ 
gible  multiple  scattering);  these  constraints  are  generally  consistent  with 
the  experimental  conditions  of  present  interest.  The  limitation,  of  course, 
will  come  in  the  assessment  of  particles  in  the  0.1  to  1  pm  size  range. 

An  alternative  technique  which  has  been  proposed  for  particle  sizing  in 
the  range  of  0.02  <  a  <  1  pm  is  a  laser-modulated-incandescence  (LMI)  tech¬ 
nique  (Ref.  31).  Briefly,  for  the  LMI  method  a  short  pulse  laser  operating  at 
a  wavelength  strongly  absorbed  by  the  particle  Uabs  <  rp)  is  used  to  irradi¬ 
ate  and  heat  a  sample  of  the  particle  field.  The  time  history  of  the  result¬ 
ing  particle  incandescence  is  then  used  to  infer  the  particle  dimensions.  The 
particle  size  affects  the  particle's  initial  thermal  equilibration  by  heat 
conduction  and  its  subsequent  cooling  by  radiation  and  heat  transfer  to  the 
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surrounding  medium.  The  advantage  of  this  technique  is  that  it  is  relatively 
easy  to  implement  and  can  be  used  for  assessing  particle  sizes  in  a  regime 
inaccessible  to  many  other  in-situ  optical  techniques. 

4.2.2  Diagnostic  Options  for  Extinction,  Emission,  and  Scattering 
Measurements 


4.2.2. 1  Radiation  Sources  for  Extinction  (and  Scattering)  Measurements 

Line  Sources  -  The  sources  considered  in  this  category  basically  consist  of  a 
variety  of  commercial  laser  sources  and  atomic  discharge  lamps.  In  general, 
lasers  are  not  needed  for  extinction  measurements  but  make  the  task  easier. 

On  the  other  hand,  lasers  are  virtually  essential  for  performing  angularly 
resolved  scattering  measurements  with  acceptable  signal-to-noise  levels. 

In  the  category  of  atomic  discharge  lamps,  one  attractive  option  is  to 
use  a  commercial  hollow  cathode  discharge  lamp  with  a  stabilized  power  supply, 
tor  example  a  llg  or  Ilg-Cd  lamp  can  provide  approximately  nine  strong  discrete 
lines  distributed  over  the  wavelength  range  of  300  to  500  nm  with  intensity 
levels  more  than  adequate  for  making  extinction  measurements  during  the  shock 
test  time. 

Commonly  available  commercial  lasers  are  also  useful  and  are  essential  it 
good  quality,  quantitative  scattering  data  is  to  be  obtained.  Attractive  can¬ 
didate  laser  sources  for  these  measurements  are  listed  in  Table  9. 

Calculations  of  the  optical  signals  that  can  be  expected  in  a  practical 
shock  tube  scattering  measurement  indicate  that  to  achieve  detectable  scat¬ 
tering  will  require  laser  powers  of  at  least  a  few  milliwatts.  This  laser 
power  requirement  increases  dramatically,  however,  if  the  particles  to  be 
probed  are  highly  emissive  and  at  high  temperature.  This  is  illustrated  in 
Figure  47  where  we  have  plotted  the  approximate  laser  power  requirements  for  a 
scattering  measurement  at  1060  nm  (A A  =  10  nm).  For  the  purpose  of  the  calcu¬ 
lation,  an  acceptable  scattering  signal  was  assumed  to  be  one  for  which  the 
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TABLE  9.  Candidate  Commercial  Lasers  Sources  for  Extinction 
and/or  Scattering  Measurements. 


Wavelength  (pm) 

Laser  Type 

0.458+0.515 

Argon -ion 

0.633 

He-Ne 

1.06 

CW  Nd : YAG 

or  Glass 

3.39 

IR  He-Ne 

9.4+10.8 

CW  co2 

signal  is  a_t  least  20  percent  of  the  background  signal  due  to  particle  emis¬ 
sion.  The  two  solid  curves  in  the  plot  represent  the  radiant  emission  that 
can  be  expected  from  the  particles  for  different  assumed  values  of  the  parti¬ 
cle  field  eraissivity  (e  =  0.01  and  e  -  0.1).  As  can  be  seen,  for  particle 
temperatures  >  1400  K  and  e  >  0.01,  the  laser  power  requirements  for  a  good 
scattering  measurement  increase  dramatically  with  the  particle  temperature. 
For  example,  for  Tp  =  2500  K  and  e  =  0.01,  the  required  laser  power  increases 
to  nearly  300  raw. 

In  general,  the  problem  of  background  particle  emission  will  be  even 
worse  at  longer  wavelengths  but  can  improve  significantly  when  working  at 
shorter  wavelengths  (such  as  at  the  468  nm  Ar-ion  line).  Figure  48  shows  the 
results  of  calculations  for  the  latter  wavelength.  Comparing  Figures  4/  and 
48,  one  can  clearly  see  the  advantage  of  operating  at  the  shorter  wavelengths. 
For  example,  assuming  a  particle  emissivity  of  0.1  at  both  wavelengths  and  a 
temperature  of  2200  K,  we  see  that  the  laser  power  required  for  a  scattering 
measurement  decreases  from  approximately  1 . 5W  at  1.06  pm  to  only  ~50  inW  at 
490  nm. 
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Figure  48.  Approximate  laser  power  requirements  for  scattering 
mea.F.9.F.e?>Pn.t  ^scatter  >  ^.2  pemission^  • 
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In  the  final  experiment  design,  similar  calculations  should  be  carried 
out  for  other  laser  wavelengths  of  interest  and  for  a  wider  range  of  particle 
emissivities.  In  this  way,  acceptable  operating  regimes  for  scattering  meas¬ 
urements  can  be  clearly  defined  prior  to  hardware  implementation  and  experi¬ 
mental  testing. 

Continuum  Sources  -  In  addition  to  using  an  array  of  discrete  line  sources  for 
these  experiments,  it  is  also  possible  to  make  use  of  broadband  continuum 
sources  for  performing  the  extinction  measurements.  For  example,  a  good 
source  for  use  in  the  wavelength  range  of  0.4  pm  <  X  <  5  pm,  would  be  a 
tungsten  ribbon  lamp  equipped  with  a  sapphire  window.  Such  sources  are 
commercially  available,  and  can  be  operated  with  filament  temperatures  up  to 
2800  K.  They  routinely  provide  spectral  radiances  that  start  from  approx¬ 
imately  5  pW/nm  sr  mm2  at  0.4  pm,  peak  to  ~150  pW/sr  mm2  nm  at  approximately 
1  pm,  and  then  fall  to  about  3  pW/sr  mm2  nm  at  4.5  pm.  Assuming  a  solid  angle 
collection  factor  of  3x10-2  sr  (12°  cone  angle)  and  an  effective  filament  area 
of  2x6  =  12  mm2,  these  radiances  yield  availble  spectral  power  densities  of 
1.7  pW/nm  at  0.4  pm,  54  pW/nm  at  1  pm,  and  1.2  pW/nm  at  4.5  pm.  In  all  cases, 
such  spectral  power  densities  should  be  more  than  adequate  for  obtaining  good 
quality  shock  tube  extinction  measurements.  Scattering  measurements,  on  the 
other  hand,  would  be  extremely  problematic  with  such  a  source. 

For  wavelengths  >4.5  pm,  a  standard  large  area  blackbody  calibration 
source  (Tmax  =  1000C)  should  be  adequate  for  obtaining  extinction  data.  For 
example,  with  T  =  1250  K  and  a  source  emissivity  =  1,  the  radiance  at  4.5  pm 
is  5.4  pW/sr  mm2  nm.  Assuming  a  source  area  =  50  mm2  and  Aflco^t  =  3x10-2  sr, 
this  yields  an  available  spectral  power  density  of  8.1  pW/nm.  If  the  particle 
loading  can  be  controlled  so  that  the  aerosol  emissivity  doesn't  exceed  0.1, 
extinction  signals  that  are  measurable  above  the  background  emission  should  be 
measurable  over  the  full  temperature  range  of  interest.  In  practice,  the  use 
of  the  1000C  blackbody  source  for  extinction  measurements  might  actually  be 
feasible  down  to  wavelengths  as  short  as  2.5  microns. 
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4. 2. 2. 2  Detection  Schemes  for  Extinction,  Scattering  and  Emission 
Measurements 


Detection  in  the  Near  UV,  Visible,  and  Near  IR  -  One  device  that  would  be 
extremely  useful  for  obtaining  spectral  information  from  the  near  UV  to  near 
IR  (\  <  1  pm)  in  a  single  shock  tube  run  is  a  ga ted-intensi f ied  diode  array 
spectrograph  (also  known  as  an  optical  multichannel  analyzer  (OMA) ) .  PSI  cur¬ 
rently  owns  such  a  device  which  can  routinely  be  used  to  obtain  spectra  over 
the  range  from  approximately  0.3  to  0.85  |jm,  with  coverage  on  a  single  shot 
presently  limited  to  a  span  of  approximately  0.4  (im*  With  0.4  j i m  coverage, 
the  spectral  resolution  is  ~20A.  The  image  intensifier  feature  provides 
excellent  sensitivity  with  overall  quantum  efficiencies  corresponding  to  a 
response  of  0.1  count-per-incident  photon  in  the  visible  and  UV.  The  gating 
option  is  also  attractive  as  it  allows  time-resolved  information  to  be 
obtained  at  selected  delay  times  after  shock  passage.  Detection  time  windows 
can  routinely  be  varied  from  as  short  as  180  ns  to  up  to  6  ms  with  control¬ 
lable  delay  times  from  270  ns  to  7  ms.  The  maximum  repetition  frequency  for 
obtaining  spectra  is  60  Hz.  Thus,  only  one  spectrum  can  be  recorded  during 
the  available  shock  test  time  of  a  few  milliseconds. 

Despite  this  single  spectrum  limitation,  the  possibility  may  still  exist 
to  record  extinction  and  emission  spectra  on  a  single  shot.  The  approach 
would  involve  using  a  combination  of  discrete  line  sources  to  measure  the 
extinction  and  superimposing  these  extinction  signals  upon  the  broadband 
emission  spectrum.  The  procedure  would  be  as  follows:  (1)  before  particle 
injection  and  shock  initiation  (and  >  16  ms  before  shock  passage)  a  spectrum 
of  the  transmitted  intensity  of  the  line  source  through  the  cold  bath  gas  is 
recorded,  (2)  then,  subsequent  to  particle  dispersal  and  shock  passage, 
another  spectrum  of  the  attenuated  line  spectrum  is  recorded  superimposed  on 
the  broadband  thermal  emission  spectrum  from  the  shock-heated  aerosol.  A  pic¬ 
torial  representation  of  the  anticipated  spectra  is  given  in  Figure  49. 

Another  approach  to  obtaining  emission,  extinction,  (and  scattering)  data 
in  the  near  UV  to  near  IR  is  to  use  an  array  of  dedicated  photomultiplier 
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Figure  49.  Anticipated  OMA  spectra  for  combined  extinction  and 
emission  measurements. 


tubes  and/or  solid  state  detectors  each  assigned  to  a  fixed  wavelength  band¬ 
pass.  With  this  detection  scheme,  simultaneous  emission  and  extinction  (or 
scattering)  data  can  be  obtained  by  high  speed  light  choppers  and  various 
laser/line  sources.  The  scheme  is  illustrated  schematically  in  Figure  50  with 
a  representation  of  the  expected  temporal  waveform  given  in  Figure  51. 

The  one  advantage  of  this  multiple  independent  detectors  approach  is  that 
it  provides  a  more  complete  time  history  of  the  aerosol  extinction  and  emis¬ 
sion  behavior  after  shock  passage.  The  disadvantage  is  that  it  requires  many 
separate  detector  and  data  recording  channels  and  provides  less  complete  spec¬ 
tral  coverage  than  the  diode  array  spectrograph  does.  The  best  solution  is 
perhaps  to  combine  one  or  two  independent  eraission/extinction  bandpass  detec¬ 
tors  with  the  diode  array  spectrograph  scheme  described  above. 

Another  spectral  region  of  interest  is,  of  course,  the  mid-IR 
(A  =>  2+5  um)»  In  this  region,  it  is  probably  best  to  use  a  few  dedicated  InSb 


71 


DETECTOR 


SOURCE 


Figure  50.  High  speed  chopping  scheme  for  recording  extinction  and  emission 
signals  on  a  single  detector. 
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Figure  51  .  Anticipated  temporal  waveform  with  high  speed  chopping  scheme 
for  combined  emission/extinction  measurements  in  shock  tube. 
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detectors  each  equipped  with  a  different  bandpass  filter.  Bandpass  centers  of 
interest  might  be,  for  example,  2.8,  3.39,  4.3  and  4.8  pm.  The  3.39  n m  band¬ 
pass  is  attractive  because  there  is  an  IR  He-Ne  laser  available  at  this  wave¬ 
length  which  can  be  used  to  simultaneously  measure  extinction  and  emission. 

It  may  also  be  possible  to  use  the  tungsten  ribbon  lamp  with  sapphire  window 
described  earlier  to  measure  extinction  at  other  wavelengths. 

The  final  wavelength  regime  of  interest  for  the  present  measurements  is 
what  we  will  refer  to  here  as  the  long  wave  IR  region,  or  \  =  0  pm  ■+■  12  pm. 

For  spectral  coverage  in  this  regime  we  recommend  using  a  multielement  linear 
array  of  liquid  nitrogen  cooled  HgCdTe  detectors  coupled  to  a  low  resolution 
IR  spectrograph.  For  example,  with  an  eight  element  array  it  should  be  pos¬ 
sible  to  obtain  a  spectrum  from  6  to  11  pm  in  a  single  shot  with  a  resolution 
of  approximately  0.6  pm.  Continuous  spectral  coverage  in  this  wavelength 
regime  appears  to  be  important  because  of  the  onset  of  IR  absorption  bands  for 
several  of  the  materials  of  interest.  Once  again,  for  obtaining  simultaneous 
extinction  and  emission  data  in  this  regime  it  should  be  possible  to  use  a 
conventional  blackbody  source  (Tmax  =  1000°C)  with  the  high  speed  chopping 
technique  described  earlier.  For  the  latter,  the  detector  channels  could  be 
either  multiplexed  or  read  out  independently.  The  latter  would  permit  a  full 
time  history  to  be  obtained  for  each  wavelength  channel. 
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5.  CONCLUSIONS 


After  considering  several  potential  techniques  for  producing  the  desired 
range  of  particle  environments,  we  have  chosen  Lhe  shock  tube  method  for  the 
temperature  range  iOOO  to  3000  K.  While  it  is  desirable  to  also  obtain  data 
between  300  and  1000  K,  this  temperature  range  is  best  achieved  in  a  slow- 
flow,  variable-temperature  tube  furnace. 

The  proposed  experiments  will  be  performed  behind  reflected  shock  waves 
in  a  4-in.  i.d.  shock  tube  facility  consisting  of  a  driver  section  4m  in 
length  and  a  modular  test  section  7.2m  in  length.  Based  on  working  times 
actually  observed  as  well  as  times  calculated  from  Kankine-Hugoniot  relations, 
we  expect  working  times  ranging  from  2.5  ms  at  3000  K  to  7.2  ms  at  1000  K  for 
a  dilute  suspension  of  particulates  in  Ar.  These  times  are  well  in  excess  of 
the  predicted  times  for  convective/conductive  heating  and  equilibration  of 
0.1  to  10  um  (radius)  particles  in  Ar  at  20  atm. 

Prior  to  initiation  of  each  shock  wave  experiment,  selected  particles 
will  be  dispersed  into  the  observation  volume  near  the  end  wall  by  means  of  a 
radial  injection  system.  This  approach  offers  several  advantages  over  other 
methods  which  have  proven  successful  in  other  "dusty"  shock  tube  experiments. 
The  injector  will  be  capable  of  dispersing  ~  O.lg  samples  of  particles 
uniformly  into  the  optical  field  of  view  in  a  reproducible  way.  To  insure 
adequate  entrainment  and  mixing  characteristics,  the  jet  conditions  will  be 
optimized  in  bench  tests  using  visiible  light  scattering  to  visualize  the 
particle  flow  field. 

The  experimentally  observed  light  extinction,  particle  emissivity,  and 
scattered  light  angular  distributions  are  directly  related  to  the  effective 
cross  sections  (or  efficiency  factors)  for  extinction,  absorption,  and  scat¬ 
tering.  These  efficiency  factors  are  in  turn  dependent  on  particle  size  and 
wavelength,  and  can  be  related  to  the  complex  refractive  index  at  a  given 
wavelength  using  Mie  scattering  theory. 
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The  sensitivities  of  the  observable  extinction,  absorption,  and  scat¬ 
tering/polarization  parameters  to  their  real  and  imaginary  components  have 
been  discussed.  from  this  analysis  we  have  shown  that  ultraviolet  extinction 
measurements  may  be  used  to  determine  particle  number  densities  multiplied  by 
cross-sectional  area  (Npur^)  with  reasonable  accuracy,  while  both  extinction 
and  emission  become  sensitive  to  refractive  index  of  size  parameter  of  1.0. 
Once  the  quantity  Npirr^  has  been  determined  to  within  some  error  bound,  then 
the  quantities  and  y^  can  be  determined  at  longer  wavelengths  to  similar 
levels  of  accuracy.  This  in  turn  allows  the  determination  of  n  and  k,  within 
appropriate  error  limits,  from  computed  sensitivity  curves. 

Absolute  scattering  measurements  are  difficult  to  implement  owing  to: 

1)  the  highly  structured  angular  dependence  for  intermediate  to  large  size 
parameter;  2)  the  need  to  calibrate  the  scattering/detection  geometry  in  situ; 
and  3)  the  requirement  for  an  intense  laser  source  at  each  wavelength  of 
interest.  However,  relative  scattering  measurements  at  two  or  more  observa¬ 
tion  angles  can  be  extremely  important  as  a  cross-check,  of  the  preshock  parti¬ 
cle  size  determinations  and  as  a  probe  for  agglomeration  effects.  Similarly, 
the  polarization  ratios  of  scattered  light  provide  a  sensitive  test  of 
particle  shape  effects. 

Finally,  the  effects  of  particle  size  distribution  and  agglomeration  have 
been  shown  to  significantly  affect  the  measured  extinction  and  emission 
efficiencies.  This  effect  mandates  careful  characterization  of  the  particle 
size  distribution  for  the  experiments. 

For  the  basic  radiometric  diagnostic  approach,  we  recommend  simultaneous 
measurements  of:  1)  light  extinction  at  several  selected  wavelengths  in  the 
ultraviolet,  visible,  and  infrared,  using  a  variety  of  light  source/  detector 
combinations;  2)  spectrally  resolved  particle  emission  in  the  visible  red  and 
i-frared  using  calibrated  radiometric  and/or  spectrometric  detection  schemes; 
and  3)  laser  scattering  angular  dissymmetry  ratios  and  polarization/  depolari¬ 
zation  ratios  tor  confirmation  of  particle  size  and  shape  characteristics. 
Laser  scattering  and  crossed-path  extinction  measurements  will  also 


be  used  to  characterize  the  spatial  distributions  of  the  dispersed  particles 
in  the  reflected  shock. 

The  recommended  program  of  experimental  measurements,  coupled  with  the 
static  rocket  firing  program  at  AFRPL,  will  lead  to  the  development  of  a 
unique  data  base  on  particulate  radiative  properties  which  is  sorely  needed 
for  predictive  modeling  of  plume  signatures.  Possible  extensions  of  this  work 
would  include  more  in-depth  studies  of  the  mechanisms  influencing  particulate 
tormation  and  contamination  in  plumes,  development  of  means  tor  "tailoring" 
the  particle  radiation  by  manipulating  these  mechanisms,  and  application  of 
our  diagnostic  methods  to  lower  temperature  measurements  (300  to  1000  K.)  in  a 
variable-temperature  flow  tube. 
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